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I The Cryogenic Dark Matter Search (CDMS-II) employs low-temperature Ge and Si detectors to 

, seek Weakly Interacting Massive Particles (WIMPs) via their elastic scattering interactions with 

J — ' nuclei. Simultaneous measurements of both ionization and phonon energy provide discrimination 

' against interactions of background particles. For recoil energies above 10 keV, events due to back- 

l/^ , ground photons are rejected with > 99.99% efficiency. Electromagnetic events very near the detector 

c/j ■ surface can mimic nuclear recoils because of reduced charge collection, but these surface events are 

^5 ' rejected with > 96% efficiency by using additional information from the phonon pulse shape. Ef- 

, ficient use of active and passive shielding, combined with the the 2090 m.w.e. overburden at the 

II ■ experimental site in the Soudan mine, makes the background from neutrons negligible for this first 

Q ' exposure. All cuts are determined in a blind manner from m situ calibrations with external radioac- 

^ \ five sources without any prior knowledge of the event distribution in the signal region. Resulting 

. efficiencies are known to ~10%. A single event with a recoil of 64 keV passes all of the cuts and is 

' consistent with the expected misidentification rate of surface-electron recoils. Under the assump- 

^ ' tions for a standard dark matter halo, these data exclude previously unexplored parameter space for 

• 1— ( , both spin-independent and spin-dependent WIMP-nucleon elastic scattering. The resulting limit on 

' the spin-independent WIMP-nucleon elastic-scattering cross-section has a minimum of 4x lO"**^ cm^ 

$H ' at a WIMP mass of 60 GeV c-^. The minimum of the limit for the spin-dependent WIMP-neutron 
elastic-scattering cross-section is 2 x 10"'^^ cm^ at a WIMP mass of 50 GeV c~^. 

PACS numbers: 95.35.-|-d, 95.30.Cq, 85.25.0j, 29.40.Wk, 14.80.Ly 



I. INTRODUCTION 



This paper reports new results from the Cryogenic 
Dark Matter Search (CDMS-II) experiment, a search for 
non-luminous, non-baryonic Weakly Interacting Massive 
Particles (WIMPs) 1] that could form the majority of 
the matter in the universe 2, 3J. The scientific case for 
WIMPs continues to grow stronger, most recently com- 
bining the WMAP results with studies of large-scale 
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clustering, the Sloan Digital Sky Survey 0,01 j and su- 
pernova redshift data 0,13 ■ 

One attractive WIMP candidate is the lightest su- 
persymmetric particle (LSP), which arises naturally in 
many models of supersymmetry |^ IT(il | and is stable if 
R-parity is conserved [ij. Some of these recent models 
of supersymmetry favor a LSP mass in the range 50- 
500 GeV c-2. 

A reasonable model for the distribution of WIMPs in 
our own galaxy that is consistent with measurements of 
spiral galaxy rotation curves 12] is a roughly isothermal 
spherical halo around our galaxy with a mean velocity of 
^ 230 km s~^. With this mass range and velocity, the 
kinetic energy imparted to a nucleus in an elastic WIMP- 
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nucleon scattering event Il3l Il4| would range from a 
few keV to tens of keV [l^. The expected event rate 
(<1 event keV~^ kg^^ day""'^) is estimated from previous 
experimental exclusion limits, inferences from the anni- 
hilation cross-section consistent with the present WIMP 
relic abundance estimate, and the latest supersymmetric 
models. 

The small recoil energy, coupled with the expected low 
event rate, means that it is vital to suppress backgrounds 
effectively. Active and passive shielding are used to re- 
duce backgrounds produced outside of the experimental 
apparatus, leaving decays of radioactive contaminants in- 
side of the shielding as the dominant natural background. 
The products of these decays all interact electromagnet- 
ically, so that it is particularly useful to discriminate 
between electron-recoil events (most backgrounds) and 
nuclear-recoil events (WIMPs and background neutrons). 

The CDMS-II experiment discriminates nuclear recoils 
from electron recoils by measuring both the ionization 
and phonon energies of interactions within Ge and Si 
detectors. The simultaneous measurements provide a 
primary discrimination of recoil type, as the ionization 
signal for nuclear recoils is suppressed relative to elec- 
tron recoils. The simultaneous measurements also al- 
low an accurate determination of the recoil energy for 
both nuclear and electron recoils. The pulse shape of the 
phonon signal gives further discrimination against events 
that could potentially be misidcntificd as nuclear recoils 
because of incomplete charge collection near the detec- 
tor surface |l6l IT^ . These events are primarily due to 
low-energy electrons. 

The remaining background comes from neutrons, 
which produce nuclear-recoil events identical to WIMPs 
and so must be distinguished by other means. First, neu- 
trons often scatter several times within the detector ar- 
ray, while a single WIMP will not. Second, while Ge and 
Si have similar scattering rates per nucleon for neutrons, 
the WIMP-nuclcon scattering rate is expected to be 5-7 
times greater in Gc than in Si for all but the lowest-mass 
WIMPs. Third, the kinematics of neutron elastic scat- 
tering give a recoil energy spectrum scaled in energy by 
a factor of ^ 2 in Si compared to Ge, whereas the factor 
would be ~1 or less for WIMP elastic scattering. The 
multiple-scatter consideration can be used to identify a 
fraction of the WIMP candidate events as being uniquely 
due to neutrons. All three of these methods can be used, 
in conjunction with Monte Carlo simulations, to statisti- 
cally subtract any neutron background. 

The CDMS collaboration is operating the CDMS-II ex- 
periment in the Soudan Underground Laboratory (Min- 
nesota, U.S.A.) [ill. The 780 m (2090 meters water 
equivalent) of rock overburden reduces the surface muon 
flux by a factor of 5 x lO''. In the shallow underground 
site at Stanford University 0,113, hereafter referred to 
as SUF, neutrons produced by cosmic-ray muons had be- 
come the limiting background of the CDMS-I and -II ex- 
periments. With this background reduced by a factor 
of 300 at Soudan, the existing CDMS detectors have 



a greatly improved sensitivity to rare WIMP-scattering 
events. 



This paper presents results from the first dark mat- 
ter run taken in the Soudan mine, from October 2003 
through January 2004, with the same six detectors pre- 
viously run at the SUF shallow site (HIHIlll. For the 
first dark matter run of these six detectors at the Soudan 
mine, we will describe in this paper the methods and re- 
sults of the two analyses used. The first analysis is 
referred to here as the "initial" analysis and the second 
is referred to as the "current" analysis. Both analyses 
use identical cuts which were set without examining the 
WIMP-search data. This blind method of setting cuts 
guarantees that the cuts are determined without bias 
from details of the WIMP-search data. The two anal- 
yses of the WIMP-search data differ slightly in the algo- 
rithm used to estimate the ionization pulse height (see 
Sec. IV AT)) . In the initial analysis, the ionization pulses 
of approximately half of the data were inadvertently an- 
alyzed using a time-domain algorithm with poorer reso- 
lution in one of the ionization channels (see Sec. IV A"I|l . 
The loss of resolution caused one of the cuts to be un- 
intentionally overly severe (see Sec. IVB4|) . leading to a 
shght reduction in our exposure (see Sec. lVCjl . This ini- 
tial analysis found no WIMP candidates in the WIMP- 
search data. The current analysis uses the intended al- 
gorithm for analyzing all the ionization pulses. Since the 
current analysis was conducted after the WIMP-search 
data had been studied using the initial analysis, we do 
not consider the results from the current analysis to be 
blind in a formal sense. We emphasize, however, that 
since the cuts for the current analysis were set before 
examining the WIMP-search data with either analyses 
described, this later analysis still remains unbiased. The 
current WIMP search analysis finds one WIMP candi- 
date event, but this is also consistent with the expected 
background. For either analysis, the lack of a signal rules 
out a significant new range of WIMP models. 



Section |n] of this paper describes the CDMS-II ex- 
perimental apparatus at Soudan, including the cryogen- 
ics, shielding, veto, ZIP detectors, readout electronics, 
and data acquisition systems. Section IIIII presents an 
overview of the experimental run, triggering, and cali- 
brations with external sources. Section IIVI reports on 
the neutron, photon, beta and alpha backgrounds ob- 
served during the WIMP-search run and comparisons 
with expectations from Monte Carlo simulations. Sec- 
tion^presents the analysis of our data. In that section, 
we describe our pulse-shape analysis, event reconstruc- 
tion, WIMP selection cut definitions and efficiencies, es- 
timates of our exposure and backgrounds, and estimates 
of our systematic uncertainties. Section IVII presents the 
results and discusses the importance of our new WIMP- 
nucleon cross-section exclusion limits. 
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II. THE SOUDAN INSTALLATION 

A. The Soudan Underground Laboratory 

The CDMS-II experiment is housed in one of the two 
excavated caverns of the Soudan Underground Labora- 
tory, operated by the University of Minnesota. The lab- 
oratory is located in the Soudan Underground Mine State 
Park, operated by the Department of Natural Resources 
(DNR) of the State of Minnesota. Access to the mine is 
provided by DNR personnel, and the laboratory techni- 
cal and management staff is employed by the University 
of Minnesota. The CDMS infrastructure in the Soudan 
Underground Laboratory was designed and constructed 
jointly by Fermilab and the University of Minnesota. The 
laboratory is situated at a depth of 780 m, and muon flux 
measurements indicate a water equivalent depth of about 
2090 m. 



B. Cryogenics 

The CDMS-II detectors operate at a temperature of 
50 mK inside a specially made cold volume cooled by 
an Oxford Instruments Kelvinox 400-S dilution refriger- 
ator. The cryogenic system has followed the same basic 
design used in the original CDMS low-background facil- 
ity at Stanford University [2(1 I21I l23 | , but with many 
additions to allow automatic and remote control. These 
improvements make it possible to maintain a stable and 
robust detector environment even when physical access 
to the laboratory is limited or unavailable. 

The dilution refrigerator, attached cold volume, and 
surrounding shielding for CDMS-II are housed in an RF- 
shielded enclosure at the Soudan Underground Labora- 
tory. The enclosure, known as the RF room, has been 
measured to be a class-10,000 clean room during work- 
ing hours, and better than a class-1000 clean room when 
unoccupied. The pumps, cryogen supplies and cryogenic 
control systems for the dilution refrigerator are situated 
outside of the RF room in the "cryo-pad" area. The de- 
tector front-end electronics are in crates inside the RF 
room, adjacent to an electrical break-out box attached 
to the detector cold volume. All other readout electron- 
ics for the detectors and veto system are housed outside 
the RF room in an electronics room. 

The detector cold volume is the innermost of six 
nested, cylindrical cans that together make up the CDMS 
cryostat or "icebox" [S^. These cans are made 
from OFHC, low-radioactivity copper and are individu- 
ally thermally coupled to the various temperature stages 
of the dilution refrigerator through a horizontal configu- 
ration of five nested copper tubes and one solid cold fin- 
ger, collectively called the "cold stem." The innermost 
can, connected to the mixing chamber of the refrigera- 
tor by the solid copper finger at the center of the cold 
stem, is 30 cm in diameter and 30 cm high, and can ac- 
commodate up to seven stacks, or "towers," of six ZIP 



detectors each at < 50 mK. The detector volume is cooled 
efficiently by conduction through the copper stems but 
kept away from all cryogenic liquids and any radioactiv- 
ity in the commercial refrigerator itself. A separate hor- 
izontal set of copper stems, called the "electronics stem" 
contains striplines connecting individual detector stacks 
to the room temperature Front-end Electronics Boards 
(FEBs) in the RF room. 

Only radioactively-screened, low-background materials 
were used in the fabrication of the icebox. The icebox is 
surrounded by a 2-mm-thick mu-metal shield at room 
temperature which screens the detectors and cold elec- 
tronics from ambient magnetic fields. 

A Moore APACS industrial control unit controls much 
of the cryogenic system. It continuously monitors tem- 
peratures, pressures, and flows in the cryogenic system, 
all of which are stored and displayed via commercial soft- 
ware (Intellution) running on a PC and are also acces- 
sible via a web interface. The APACS unit has been 
programmed via the Intellution software to automati- 
cally perform liquid helium and nitrogen transfers, and 
to respond to certain situations that could otherwise re- 
quire operator intervention. Whenever any of the mea- 
sured temperatures, pressures, liquid levels, or flow rates 
fall outside normal bounds, the operator is contacted by 
phone and can remotely resolve most problems using a 
secure computer interface to APACS. A separate sys- 
tem, the Oxford Instruments "Intelligent Gas Handling" 
(IGH) unit, monitors and directs the circulation of the 
■^He/'^He mixture. This system, which uses LabVIEW 
software, also supports remote monitoring and operation. 
Uninterruptible Power Supplies (UPSs) provide stability 
when electrical power is interrupted, as happens a few 
times a year at Soudan during summer thunderstorms. 

Except for a few brief cryogenic excursions (described 
in ) , the six detectors described here were kept at base 
temperature {r^ 50 mK) from July of 2003 until August, 
2004. During this entire run, the dilution refrigerator 
suffered from a small leak from the liquid helium reser- 
voir into the insulating (outer) vacuum space. Two tur- 
bomolecular pumps continuously evacuated this vacuum 
space, thereby maintaining an acceptably low pressure 
for thermal isolation. Nonetheless, helium constantly ac- 
cumulated on the surfaces at 4 K inside the refrigerator 
and the icebox. Eventually (about once every month), 
enough He accumulated to cause a sudden softening of 
the vacuum (to about 5 x 10~^ torr). This caused a sub- 
stantial boil-off of liquid helium in a very short time pe- 
riod (<1 min.), and temporarily raised the temperature 
of the detectors as high as 1 K. The leak also increased 
the average rate of helium boil-off so that two cryogen 
transfers were required each day. Detector operation was 
suspended during these transfers, reducing our detector 
operation live time fraction. The leak did not expose the 
detectors to helium, because they were inside a separate, 
inner vacuum space. 

The icebox also contains the cold hardware needed to 
mount and operate the detectors. Each detector tower 
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has four temperature stages (4 K, 600 mK, 50 mK, base) 
and contains the wiring and heat-sinking required for six 
detectors, which are suspended below and are stacked 
with 2 mm separation between neighboring detectors. 
The six detectors used in this WIMP-search together 
with their tower are collectively referred to as "Tower I." 
Coaxial wires passing through the tower connect each 
detector to its own cold electronics card mounted on 
top of the tower. This card contains two Field Effect 
Transistors (FETs) for the readout of the two ioniza- 
tion channels. The FETs are weakly heat-sunk to the 
4-K stage, and self-heat to 130 K for optimal noise per- 
formance. A separate card at 600 mK contains the four 
arrays of Superconducting Quantum Interference Devices 
(SQUIDs) required for the readout of each detector's four 
phonon channels. The signals are carried to room tem- 
perature through flexible, 2.5-cm-wide, 3-m-long copper- 
kapton striplincs. These are heat-sunk at 4 K and 77 K, 
and pass through a copper radiation shield as they pro- 
ceed through the electronics stem into a breakout box 
at room temperature. All these fixtures have been de- 
signed to minimize the heat conducted and radiated into 
the cold stages, and to allow the FETs to operate with 
suitably low noise. The cold hardware is constructed 
from radioactively-screened, low-background materials, 
in particular copper, kapton, and custom low-activity sol- 
der [23. 



C. Shielding and Muon Veto 

To minimize the number of interactions that could im- 
itate nuclear recoils and to decrease the overall back- 
ground rate, the CDMS-II experiment employs several 
layers of passive shielding as well as an active muon veto. 
Backgrounds that need to be shielded against include 
high-energy neutrons produced by cosmic ray muons in- 
teracting with the surrounding rock or shield; gamma- 
rays and neutrons from radioactivity in the surrounding 
rock; and photons, electrons, and alphas from radioactive 
impurities on surfaces within the lead shielding. 

The CDMS-II experiment at the Soudan Underground 
Laboratory is situated beneath an overburden of 2090 
meters water equivalent (m.w.e.). This overburden re- 
duces the surface muon flux by a factor of 5 x 10*. Within 
the RF room, the layers of shielding are arranged con- 
centrically around the icebox (see Fig. Outermost is 
an active muon veto which tags the remaining muon flux. 

Forty scintillator panels surrounding the passive shield- 
ing and icebox compose the active muon veto system. 
Each panel consists of a 5-cm-thick slab of Bicron BC- 
408 plastic connected to one or two 2-inch Hamamatsu 
R329-02 photomultiplier tubes. The veto panels are ar- 
ranged so that adjacent panels overlap. The top panels 
extend well beyond the side panels, guaranteeing that the 
small gap between top and side panels does not allow a 
direct line of sight to the detectors inside the icebox. 

The veto system can distinguish between muons and 



ambient photons, since minimum ionizing muons typi- 
cally deposit 10 MeV while the highest energy photon 
from natural radioactivity present is 2.6 MeV. The de- 
tection efficiency of a veto paddle approaches 100% for > 
2 MeV energy deposition. The in situ measured muon de- 
tection efficiency for the entire veto system is 99.4 ±0.2% 
for stopped muons and 99.98 ± 0.02% for through-going 
muons, where the low muon rate limits the precision of 
this measurement. On average, one muon per minute is 
incident on the veto, and the combined veto rate (domi- 
nated by ambient gammas) is ~ 600 Hz. Events with de- 
tector signals are recorded regardless of whether or not 
there was veto activity, but a subsequent analysis cut 
rejects events whose veto activity preceded the detector 
signal by less than 50 /is. A veto trigger on events for 
which multiple scintillator panels have coincident activity 
is highly efficient at identifying incident muons. 

The veto calibration is checked twice each day. Since 
there is no significant rate of muons, we calibrate the veto 
using light from a bank of blue Light Emitting Diodes 
(LEDs) fed via optical fibers to each veto panel. The 
light from the LED bank is tuned to match the scintilla- 
tion light created in each panel from a muon. Addition- 
ally, we continually monitor the veto trigger rate and the 
voltage and current supplied to each veto panel's photo- 
multiplier tube. We found no significant variation in the 
performance of the veto system over the course of the 
run. 

As depicted in Fig. ^ within the volume enclosed by 
the muon veto panels, a 40 cm thick cylindrical outer 
polyethylene layer first moderates low-energy neutrons 
from radioactive decays to below-threshold energies. In- 
side this polyethylene lies a 22.5 cm thick cylindrical lead 
shield, of which the inner 4.5 cm thickness consists of 
ancient lead [l^. The lead shielding is constructed in a 
modular way, with gaps that do not align with the detec- 
tors. An inner 10 cm thick cylindrical polyethylene layer 
provides further neutron moderation. The icebox cans 
and cold hardware provide an average shielding thick- 
ness of about 3 cm of copper directly surrounding the 
detectors. The shield is 99% hermetic, with the only 
penetrations being the cold stem and the electrical stem. 

Beginning on November 11, 2003 the air volume be- 
tween the outermost icebox copper can and the mu-metal 
shield has been continuously purged to dramatically re- 
duce the concentration of environmental radon in the 
vicinity of the detectors. The purge gas used is medical 
grade breathing air that has been stored in metal cylin- 
ders for at least two weeks, to allow most of the radon 
gas to decay. 



D. CDMS ZIP Detectors 

Each CDMS ZIP (Z-dependent Ionization- and 
Phonon-mediated) detector is a cylindrical high-purity 
Ge or Si crystal that is 1 cm thick and 7.6 cm in diam- 
eter. A single Ge (Si) ZIP has a mass of 250 g (100 g). 




SctnlitlnUir 




FIG. 1: Top view and side view of the CDMS-II shielding and veto. The detector volume is referred to as the "icebox." As 
shown, the stem to the right of the detector volume is the "cold stem" and connects the detectors and the copper cans to the 
cryostat. The stem to the left of the detector volume is the "electronics stem" and contains the wiring that connects the cold 
electronics to the room temperature electronics. 



Two concentric ionization electrodes and four indepen- 
dent phonon sensors are photolithographically patterned 
onto each crystal. 

The data described here were obtained with a single 
tower of six ZIP detectors. Within a tower, the six ZIP 
detectors are stacked 2 mm apart with no intervening 
material. This close packing not only shields the detec- 
tors from low-energy electron sources on surrounding sur- 
faces but also increases the probability that a background 
event in one detector would multiply-scatter into another 
detector. Division of the electrodes into an annular outer 
"guard" electrode and a disk-shaped inner electrode de- 
fines an inner fiducial region that is further shielded from 
low-energy electron sources or x-ray fluorescence. 

An external particle scattering in a ZIP detector can 
interact with an electron (or electrons) in the crystal {e.g. 
by Compton scattering, K-capture, etc., any of which we 
say causes an "electron recoil"), or with a nucleus (called 
a "nuclear recoil"). The interaction deposits energy 
into the crystal through charge excitations (electron-hole 
pairs) and lattice vibrations (phonons). A CDMS ZIP 
detector measures both the ionization and the phonon 
energy for every event. The simultaneous ionization and 
phonon measurement not only allows an accurate mea- 
surement of the recoil energy independent of recoil type. 



but also distinguishes between these two types of recoils. 

Depending on the material and the type of recoil, 6% - 
33% of the recoil energy is first converted into ionization 
before subsequent conversion to phonons. On average, 
one electron-hole pair is produced for every e « 3.0 eV 
(3.8 eV) of energy deposited by an electron recoil in Ge 
(in Si). The "ionization energy," Eq, is defined for con- 
venience as the recoil energy inferred from the detected 
number of charge pairs, iVq, by assuming that the event is 
an electron recoil with 100% charge-collection efficiency: 



Ec 



Q 



(1) 



Ionization energy is usually reported in units of "keVee," 
keV of the equivalent electron recoil determined from 
electron-recoil calibration measurements. The dimen- 
sionless ionization yield parameter, y, is the ratio of ion- 
ization energy to true recoil energy, y = Eq/E^, where 
is the full recoil energy for the event. This defini- 
tion of yield gives unity for electron-recoil events with 
complete charge collection measurement. 

Nuclear recoils produce fewer charge pairs, and hence 
less ionization energy, Eq , than do electron recoils of the 
same recoil energy (see Figs. |21and|31). The ionization 
yield y for nuclear-recoil events depends on both the ma- 
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FIG. 2: Ionization energy versus recoil energy for detector 
Z5 (Ge). Black dots correspond to calibration events from a 
^^^Ba source (emits gammas only) and gray dots correspond 
to calibration events from a ^^^Cf source (emits gammas and 
neutrons) . 

terial and the recoil energy, with y ~ 0.3 [y ^ 0.25) in Ge 
(in Si) for Eji > 20 keV. The simultaneous measurement 
of ionization and recoil energy therefore makes it possi- 
ble to identify and reject most electron-recoil background 
events. 



1. The Ionization Measurement 

On each ZIP detector, metal electrodes on the two 
faces of the crystal substrate serve as the sensors for the 
ionization measurement. As shown in Fig.^ the charge 
electrodes are divided into an inner fiducial electrode and 
an outer guard ring. The disk-shaped inner electrode has 
a diameter of 69 mm and so covers 82-85% of the ioniza- 
tion side of the detector. The annular outer guard ring 
is 2.0-2.7 mm wide, and a 1-mm-wide annular gap lies 
between the two electrodes. The guard ring is used to 
reject events occurring near the bare unpolished edges of 
the crystal. We reject events near these edges for three 
reasons: the ionization signal may be degraded since the 
electric field is not uniform in this region, the phonon 
response is worse since this region is not well covered 
by the phonon sensors, and background interactions are 
more likely due to the absence of self-shielding by other 
detectors in the stack. 

The electrodes are used to apply an electric field of 
a few volts/cm across the crystal, which drifts elec- 
trons and holes produced by an interaction to the de- 
tector faces. The applied electric field needs to be rela- 
tively low in order not to contaminate the phonon sig- 
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FIG. 3: Ionization yield versus recoil for detector Z5 (Ge). 
Black dots correspond to calibration events from a ^^^Ba 
source (emits gammas only) and gray dots correspond to cal- 
ibration events from a ^^^Cf source (emits gammas and neu- 
trons). The upper distribution of events are bulk electron re- 
coils which define the "electron-recoil band." The lower distri- 
bution of events are nuclear recoils which define the "nuclear- 
recoil band." 

nal with too many Neganov-Trofimov-Luke phonons (see 
Sec. Ill 0*21 ■ which would compromise the detector's abil- 
ity to discriminate between electron-recoil and nuclear- 
recoil events. The number of charges collected on the 
electrodes, TVq, is proportional to the number of electron- 
hole pairs produced by the interaction weighted by the 
drift distance of each of the charges 0, |29|| . In the case 
where all of the electrons and holes drift across the crys- 
tal, A^Q (see Eq. ^) equals the total ionization. The 
prompt recombination of the drifted electrons and holes 
in the electrodes releases all of the recoil energy from the 
electron system into the phonon system. 

There are two cases for which the measured charge, 
A'q, underestimates the total ionization. In both of these 
cases, the charges fail to drift across the entire crystal. 
The first case is due to poor charge-space neutralization 
within the crystals. At the low operating temperatures 
required for the phonon sensors, impurity sites can be left 
with a net charge. These charged impurity sites can trap 
the drifting electrons (or holes) from a recoil event. We 
neutralize the crystals routinely by flashing LEDs, whose 
photons generate excess electron-hole pairs in the crys- 
tal. These excitations can then neutralize ionized impu- 
rity sites, reducing their trapping cross-section by several 
orders of magnitude. The Si crystals contain more impu- 
rities than the Ge, and the LED photon energy spectrum 
is optimized for neutralizing Ge, thus neutralizing the Si 
detectors takes more flashing cycles than for Ge. For the 
data run described here, the LED flashing regime was 
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FIG. 4: Diagram of a CDMS ZIP detector, a) The "ionization side" of the detector with a large inner electrode and an outer 
guard ring electrode, b) The "phonon side," divided into four quadrants labeled A, B, C, and D, each consisting of 37 dies of 28 
QETs. The convention for the x-y axes is shown. The area outside the cells consists of a passive Al/W grid that is patterned 
sparsely (10% area coverage) to minimize athermal-phonon absorption while maintaining field uniformity for the ionization 
measurement, c) One of the 37 dies constituting a single phonon channel; each die contains 28 QETs. d) One of the QETs 
consisting of a l-/im-wide tungsten strip connected to 8 aluminum fins. 



not always sufficient for the Si detector Z6, which suf- 
fered incomplete charge collection for a sizable fraction of 
the data taken. The second case of underestimating the 
ionization corresponds to events occurring close to one of 
the charge electrodes. For these events, the relatively low 
applied electric drift field and self-screening from the ini- 
tial electron-hole cloud enables some of the electrons or 
holes to drift into the "incorrect" electrode. By deposit- 
ing a thin (~40 nm) layer of lightly doped amorphous sil- 
icon between each electrode and the detector surface, this 
back-diffusion is substantially reduced 0, [s^ ■ Interac- 
tions within the 10 fim "dead layer" still have deficient 
charge collection and are referred to as "surface events." 
As discussed in the following section, information from 
the athermal phonon measurement makes it possible to 
identify and reject these surface events. 



2. The Phonon Measurement 

A total of 4144 Quasiparticle-assisted Electrothermal- 
feedback Transition-edge sensors (QETs) [sj photolitho- 
graphically patterned onto one of the crystal faces form 
the phonon sensors for a CDMS ZIP. The QETs are di- 
vided into four independent channels, each consisting of 
1,036 QETs operated in parallel. Each QET consists of a 



l-/im-wide strip of tungsten (35 nm thick) connected to 
eight superconducting aluminum collection fins (300 nm 
thick), each roughly 380 /im x 55 /im, as shown in Fig.^ 

The narrow tungsten strips form the Transition-Edge 
Sensors (TESs), which are voltage biased, with the 
current through them monitored by a high-bandwidth 
SQUID array [33|- The tungsten is maintained stably 
within its superconducting to normal transition by elec- 
trothermal feedback based on Joule self-heating: if the 
sensor were hotter, the resistance would increase, de- 
creasing the current and the Joule heating; an analo- 
gous argument applies if the sensor were cooler. To en- 
sure operation in the extreme feedback limit, the sub- 
strate is kept much colder (T < 50 mK) than the 
superconducting-to- normal transition temperature of the 
tungsten sensors (Tc ~ 80 mK). Energy deposited in the 
tungsten electron system raises the temperature of the 
film, increasing its resistance and reducing the current. 

Most phonons in the crystal that reach the aluminum 
fins on the surface scatter into the aluminum, creating 
quasiparticle excitations. These quasiparticles diffuse 
through the aluminum and later enter the tungsten TES. 
Electron-electron interactions between the quasiparticles 
and the tungsten conduction electrons |33j cause an in- 
crease in the temperature of the tungsten electron sys- 
tem. Strong electrothermal feedback guarantees that the 
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power deposited into the TES is exactly compensated for 
by a reduction in Joule heating. Thus, for a voltage bias 
across the TES, Vtes, and a TES current decrement A/, 
the total energy deposited, Etes, is just 

Etes = Vtes J Aldt. (2) 

The energy in the phonon system includes not only 
the full recoil energy, Er, but also energy from drift- 
ing the electrons and holes across the crystal. The work 
done in drifting the electrons and holes manifests itself as 
Neganov-Trofimov-Luke phonons 0, . The Neganov- 
Trofimov-Luke phonons contribute to the total phonon 
energy, Ep, yielding 

Ep^En + eJ^Edq, (3) 

where the sum is taken over all charges, -Er is the recoil 
energy, E is the electric field, and dq is the distance that 
each charge drifts. Since the electric field is constant 
throughout most of the detector volume, we can rewrite 
Eq. ||2Jl as 

Ep^E^ + eV^y^^, (4) 

^do 

where do is the thickness of the crystal and Vb = Edo is 
the bias voltage between the electrodes. The summation 
in Eq. Q is the number of charges weighted by their 
drift distances and is equal to the charge collected on the 
electrodes, Nq. The definition of ionization energy given 
by Eq. (Q) gives 

Ep^En + eV^NQ ^E^ + ^Eq, (5) 

e 

where Vb is the bias voltage across the detector. Equa- 
tion (|3J) is valid even for events with incomplete charge 
collection (due, for example, to trapping or recombina- 
tion in the wrong electrode). Since we calibrate elec- 
tron recoils with full charge collection to have Eq = Ep , 
Ep = (1 + eVb/e) Eji for these events. Calibration of the 
detectors at several bias voltages using photon sources 
confirms that e w 3 eV (3.8 eV) in Ge (in Si). For elec- 
tron recoils with full charge collection in Ge at 3 V bias 
(the bias voltage for most of the data described here), 
Ep = 2Ep_. In practice, the recoil energy Ep of an event 
is inferred from measurements of the phonon and ioniza- 
tion energies: 

Ep, = Ep^ ^Eq. (6) 

Equation © is valid for all events independent of recoil 
type and the ionization collection efficiency. 

One of the advantages of measuring the athermal 
phonon signal is that it provides sensitivity to phonon 
physics that is dependent on the nature of the interac- 
tion and the event location. This sensitivity makes it 



worthwhile to discuss the phonon dynamics within the 
crystal. 

An interaction in the crystal produces high-frequency 
phonons which propagate quasi-diffusively through the 
crystal l37| . Quasi-diffusive propagation arises from 
the combination of two scattering processes: elastic scat- 
tering, which mixes phonon modes, and anharmonic de- 
cay, which increases the total number of phonons while 
reducing the average phonon frequency. The net effect 
of the quasi-diffuse propagation is a phonon "ball" ex- 
panding at roughly 1/3 the speed of sound. When the 
frequency of a phonon is sufficiently low (< 1 THz), the 
mean free path becomes comparable to the size of the de- 
tector. Such phonons are termed 'ballistic' and travel at 
the speed of sound (5 mm /xs~^ for Ge and 2.5 mm ^s~^ 
for Si) through the crystal. They rarely scatter within 
the bulk, but reflect diffusively at the crystal surfaces, 
unless they are absorbed by metal films on the surface. 
Eventually the quasi-diffuse propagation converts all of 
the initial high-frequency phonons into ballistic phonons, 
which then eventually thermalize. 

In the absence of an electric field through the crys- 
tal, no discernible difference in the phonon signal com- 
ing from electron recoils versus nuclear recoils has yet 
been demonstrated j3S |. However, during normal de- 
tector operation an electric field is present. In the first 
few hundred nanoseconds after an interaction, the freed 
electron-hole pairs drift across the entire crystal, shed- 
ding ballistic Neganov-Trofimov-Luke phonons and pro- 
ducing additional ballistic phonons upon relaxation at 
the electrodes 41]. Since an electron recoil pro- 

duces more charges than a nuclear recoil, these processes 
lead to a larger initial population of ballistic phonons 
for electron recoils. Ballistic phonons propagate at the 
speed of sound, while the energy flux of high-frequency 
phonons moves at approximately one third the speed of 
sound. This difference in the speed of propagation leads 
to a faster phonon leading edge for electron recoils when 
compared to nuclear recoils because of the larger ballistic 
fraction. 

There is an additional increase in the fraction of bal- 
listic phonons arising from the down-conversion of the 
high-frequency phonons from interactions with metal on 
the crystal surface When electrons and holes relax 
to the Fermi surface in the metal electrodes, most of the 
released energy is rapidly down-converted to a largely 
ballistic population of phonons. This effect becomes es- 
pecially significant for the special case where events are 
very close (< 1 mm) to the detector surface. For these 
events, the expanding phonon ball is sufficiently close to 
the metal on the surface that a substantial fraction of 
the initial high-frequency phonons are down-converted, 
leading to a third population of prompt ballistic phonons 
while reducing the population of high-frequency phonons. 
All of these processes produce ballistic phonons much 
more rapidly than the down conversion from quasi-diffuse 
propagation. For such surface events, the phonon signal 
will rise even faster than for bulk electron recoils. Fig- 
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FIG. 5: Peak phonon delay versus ionization yield for de- 
tector Z5. Events from a ^^^Cf source (emits both gammas 
and neutrons) are shown as gray dots. Events from a '^^^Ba 
source (emits gammas only) are shown as black dots. Surface- 
electron recoils from the ^'^^Ba calibration appear as a low- 
yield tail and having smaller peak phonon delays. Lines for 
a typical cut excluding events with high yield or low peak 
phonon delay are shown. 

ure El illustrates these effects in data. 

The measurement of athermal phonons makes it pos- 
sible to reject events in the dead layer by analysis of 
the phonon pulse shape. Additionally, it is possible to 
reconstruct the interaction location by combining the in- 
formation from the four independent phonon channels. 



E. Data Acquisition Electronics 

Custom-made electronics, the Front-end Electronics 
Boards (FEBs), are located in 9U crates near the vac- 
uum bulkhead where the striplines from the cold elec- 
tronics terminate. The back-plane of these crates define 
the "star-ground" for the entire electronics chain and is 
carefully designed to minimize noise pickup. The FEBs 
control detector settings, such as bias voltages for both 
the ionization and phonon measurements. They also con- 
tain LED flashing driver electronics, and amplifiers for 
the readout of the cold electronics' FETs and SQUIDs. 
The FEBs are programmed by GPIB commands from the 
Data Acquisition system, as indicated in Fig. 

The data acquisition (DAQ) electronics, shown 
schematically in Fig. |H| records three streams of data. 
Fast waveform digitizers (ADCs) sample the analog sig- 
nals and provide time-trace information for both ZIP de- 
tectors and veto paddles. In addition, digital time his- 
tories are recorded for veto paddle and detector chan- 



nel triggers. Finally monitoring devices (slow-sampling 
ADCs and scalers) record trigger thresholds, detector sig- 
nal offsets, and triggering rates. 

Struck SIS 3301 digitizers are used to record the ZIP 
detector signals. These sample natively at 80 MHz, but 
we average every 64 samples to achieve a sampling rate of 
1.25 MHz, which exceeds the signal bandwidth of inter- 
est. This sampling averaging gives a measured accuracy 
of 14 bits. We take 2048 samples for each of the 6 chan- 
nels (2 charge and 4 phonon) from each ZIP detector. For 
the faster scintillation pulses from the muon veto panels, 
a custom-built circuit stretches the pulses from the PMTs 
to ~ /iS widths. The signals are then waveform-sampled 
with Joerger VTR812 digitizers, operating at 5 MHz and 
taking 1024 samples. 

The trigger information for the detector analog signals 
is determined by custom-built Receiver- Trigger-Filter 
(RTF) boards, which filter the detector analog signals, 
sum the four phonon and sum the two charge channels, 
and compare the results with user-set thresholds. The 
digital output from the RTF boards is passed to custom- 
built Trigger-Conditioner Boards (TCBs), which condi- 
tion the signal. LeCroy discriminators compare the veto 
pulses with user-set thresholds (typically below 1 MeV) 
to generate the digital history for the veto system. Struck 
SIS 3400 time history buffer modules store the digital 
trigger information for both the detector and veto sys- 
tems with 1 microsecond resolution. 

A custom-built Trigger Logic Board (TLB), which pro- 
cesses the digital trigger signals for both the detector and 
veto system, determines whether the signal satisfies the 
conditions for a global trigger. There are two possible 
conditions for which the TLB issues a global trigger. The 
first condition is when there is a signal in more than one 
veto panel. This condition is designed to issue a trigger 
for muons passing through the veto. The second condi- 
tion is when the summed phonon signal of any one detec- 
tor exceeds a user-set threshold. This condition triggers 
on events in the ZIP detectors. Since the threshold for 
this criterion is fixed throughout the experiment, chang- 
ing phonon pulse-shapes (see Sec. lV A'3|l and amplitudes 
lead to variations in the location of the global trigger in 
the analog time-stream. For example, events with larger 
amplitudes (higher energy) have global triggers that are 
earlier than events with smaller amplitudes (lower en- 
ergy). Likewise, faster rising events have earlier global 
triggers than slower ones. This delay between the inter- 
action time and the global trigger can extend up to tens 
of microseconds. 



F. Data Acquisition Software 

The CDMS-II data acquisition software provides con- 
trol, monitoring, and data recording for the experiment. 
It consists of a number of independent servers written 
in JAVA and C+-t- which communicate via the JAVA 
Remote Method Invocation (RMI) and CORBA tech- 
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FIG. 6: CDMS-II data acquisition layout showing the various hardware components in relation to the detector and veto 
systems. When a detector triggers, various detector and veto elements are read out to form an event. Monitoring processes 
run continuously to monitor the health of the detector, veto, and cryogenics components. 



nologies. The use of these technologies provides a scal- 
able system with inherent remote accessibility. There 
are three important features of the software: high band- 
width, automation, and remote interfacing. Figure El de- 
picts this software system and its interface with the DAQ 
electronics. 

Though the expected event rate is low at the Soudan 
mine, the CDMS-II data acquisition software is capa- 
ble of processing data at a high event rate. Whenever 
a global trigger is generated, crate-based servers read 
out each of the VME crates that contain the digitizers 
and history buffers. This information is combined and 
recorded to disk, after which the digitizers and Trigger 
Logic Board are rearmed. The DAQ dead time is about 
20 ms, allowing us to take data at 50 Hz. This bandwidth 
provides the means to acquire calibration data at rates 



much higher than our WIMP-search data (0.1 Hz). 

During standard running most routine operations, 
such as configuration of the detector and veto electron- 
ics, periodic detector neutralization, fast veto calibration 
with the LED pulser system, data copying, and data 
archiving, are automatically performed by the data ac- 
quisition system. Additionally, if the experimental con- 
figuration deviates from the specified parameters, the 
system automatically pauses and reconfigures the exper- 
iment. This automation minimizes user intervention and 
maintains a stable experimental configuration through- 
out the entire WIMP search. 

The use of JAVA makes remote control and moni- 
toring an intrinsic part of the data acquisition system. 
This feature is useful since physical access to the un- 
derground laboratory is sometimes limited. By using 
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a thin platform-independent client, any user can moni- 
tor the CDMS-II experiment regardless of their location. 
The system broadcasts information regarding the trigger 
rates, cryogenics, high voltage, and experimental config- 
uration through data channels so that remote clients can 
access them. In addition to the remote monitoring of the 
experiment, a fast analysis of the output data, based on 
the analog and digital information from the detector and 
veto systems, is performed to produce various diagnos- 
tic plots, which are constantly updated and available on 
the web. Near real-time diagnostic plots from the on- 
line data analysis are available for review of data quality, 
including some automated features alerting operators to 
unusual performance of the apparatus. 



III. THE FIRST SOUDAN DATA RUN 
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A. Run Overview 
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The first Soudan WIMP-search run started on October 
11, 2003 and lasted until January 11, 2004. As shown in 
Fig. [7| we acquired 52.6 live-days of WIMP-search data 
with Tower 1 detectors: Zl (Ge), Z2 (Ge), Z3 (Ge), Z4 
(Si), Z5 (Ge), and Z6 (Si). The same tower, with identical 
detector configuration, had been used in a WIMP-search 
run [23, 121I l2a | at the shallow site at Stanford Univer- 
sity. The Soudan run was interrupted by brief neutron 
calibrations with a ^^^Cf source on November 25, 2003, 
December 19, 2003, and January 5, 2004. In addition, 
several extensive gamma calibrations with '^'^^Ba sources 
were performed. Most of them were taken in the first 
half of December 2003, with occasional earlier and later 
calibration runs for checks on stability. 

Daily maintenance of the experiment, including cryo- 
gen transfers into the cryostat, LED flashing periods 
that ensure detector neutralization, etc., were performed 
roughly every 12 hours, resulting in approximately 2-3 
hours per day of scheduled down-time. Occasional dif- 
ficulties with cryogenics, electronics, or the noise envi- 
ronment caused additional periods of down-time, as in- 
dicated in Fig. Hand detailed in f26| . 

B. Triggering 

Global triggers for the experiment were issued when- 
ever the sum of the four phonon traces of any detector 
exceeded a certain threshold. As described in Sec. Ill El 
the TLB determines whether to issue a global trigger by 
analyzing the digital output of all of the RTF boards and 
the veto system. This output is recorded in the digital 
history buffer. To determine the trigger efficiency for a 
given detector, we must first preselect events by requir- 
ing that the global trigger not be issued by that detector. 
The detector's trigger efficiency is the fraction of these 
preselected events for which the history buffer records 
that the detector's RTF board finds the phonon signal 



FIG. 7: Live-time before (gray) and after (black) application 
of the data-quality cut described in Sec. lV BTI The short re- 
gions of down-time correspond to episodes of electronics prob- 
lems and unintended detector warming above 1 K. 

to be above threshold within 100 /xs after the global trig- 
ger. We calculate the trigger efficiency using all of the 
WIMP-search data. For most of the detectors the trigger 
efficiency reaches 100% between 5-10 keV true recoil en- 
ergy, and for some detectors (such as Z3 in Fig. [S)) even 
lower, 3-5 keV. To achieve these efficiencies, the trigger 
thresholds were set low, such that ~ 30% of the triggers 
were caused by noise. 

C. Calibrations with External Sources 

We took calibration data on several occasions through- 
out the WIMP search using external radioactive sources. 
A ^'^^Ba source served to characterize the detector re- 
sponse to electron recoils. To avoid the Pb in the passive 
shield, the source was inserted through a special tube 
along the electronics stem of the icebox (see Fig. 
which penetrates the passive shield. This particular 
isotope offers several distinct lines: 276 keV, 303 keV, 
356 keV, and 384 keV. These lines are sufficiently en- 
ergetic that the photons can penetrate the copper cans 
of the cryostat and reach the detectors. However, only 
the ionization channels have linear response in this high- 
energy region. The phonon channels typically become 
significantly non-linear above 200 keV. Hence, we use 
these lines to calibrate the ionization channels, and then 
calibrate the phonon channels against the ionization for 
electron recoils. 

Detailed Monte Garlo simulation of the ^'^^Ba calibra- 
tion was performed using GEANT 3. The results of 
the simulation were convolved with the energy-dependent 
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FIG. 8: Phonon trigger efficiencies for some of the ZIP de- 
tectors. Black dots with error bars show the trigger efficiency 
for the detector Z3 (Ge). The other points are plotted with- 
out error bars to avoid confusion. Triangles correspond to the 
trigger efficiency for Z5 (Ge) , and diamonds correspond to the 
efficiency of Z4 (Si) which has the worst efficiency of all the 
Tower I ZIPs. 



resolution of the charge channels. The resulting predic- 
tion for Tower 1 is compared to actual calibration data in 
Fig. 1^1 The lines are very clear in all Ge detectors, mak- 
ing the calibration straightforward. For the Si detectors, 
Compton scattering dominates at these energies, making 
the lines considerably less defined. Hence, for the case of 
Si detectors, the calibration is achieved by matching the 
spectral shapes of the Monte Carlo simulation and the 
data. 

We used a ^^^Cf source to characterize the detector 
response to nuclear recoils. Figure ^1 shows a com- 
parison of the measured recoil spectrum for neutrons 
with the spectrum predicted by a simulation of the 
source. The excellent agreement demonstrates that the 
calibration of the phonon channels using electron recoils 
remains valid for nuclear recoils. This suggests that 
the phonon measurement of the recoil and associated 
Neganov-Trofimov-Luke phonons is independent of recoil 
type, or unquenched. 



IV. SIMULATIONS AND MEASUREMENTS OF 
EXPECTED BACKGROUND SOURCES 

Four ordinary particles comprise the background for 
the experiment: neutrons, gammas, betas, and alphas. 

Neutrons produce nuclear recoils and cannot be re- 
jected on an event-by-event basis unless they scatter in 
more than one detector. We simulated the expected neu- 
tron background at Soudan 43] and found that it is in- 
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FIG. 9: Comparison of the measured Ba charge spectrum 
(dots with error bars) and the Monte Carlo simulation (gray 
line) for the Ge detector Z3 (top) and the Si detector Z4 
(bottom). The expected lines for this source are at 276 keV, 
303 keV, 356 keV, and 384 keV. 



significant over the course of the run. 

Gammas and betas can be misidentified as nuclear re- 
coils when they scatter in the dead layer near the detector 
surface. The poor charge collection in the dead layer can 
lead to ionization yields that are low enough to "leak" 
into the nuclear- recoil bands, as shown in Fig. O These 
surface events (also called leakage events) comprise the 
most significant background for the experiment. The ma- 
jority of the gamma background interacts in the bulk of 
the detector where charge collection is complete. Com- 
parison of the WIMP-search bulk-electron recoils to sim- 
ulation spectra allows us to identify the sources of our 
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FIG. 10: Comparison of measured ^^^Cf neutron recoil spec- 
trum (dots with error bars) and Monte Carlo simulation (gray 
line) for coadded Ge detectors (top) and Si detectors (bot- 
tom) . The excellent agreement suggests that the phonon mea- 
surement of the recoil is unquenched. 



background gammas. 

Betas arising from contamination on the detector sur- 
face are the most difficult background to characterize. 
There are a number of possible beta emitters such as 
4°K, "c, and 2iopb. The detector Z6 has a known "c 
contamination and was included in Tower 1 as an active 
veto. We have simulated the depth distribution of vari- 
ous sources of surface events. Combining the simulation 
results with a model of the detector response shows that, 
for our analysis, betas from surface contamination are 
the dominant electromagnetic background. 



In addition to emitting betas, contaminants on the de- 
tector surface can also emit alphas. Alphas have large 
recoil energies (>1 MeV) and low ionization energies (< 
1 MeVee), making them easy to identify (see Fig. Ill|l . 
The recoiling nucleus of a typical alpha decay in our de- 
tector produces a nuclear recoil with energy of order 20- 
100 keV. Events produced by the recoiling nuclei of alpha 
decays can be mistaken for WIMPs if the emitted alpha 
particles are not detected. Our detection efficiency for 
alphas is high, so we expect this background to be neg- 
ligible. Detecting alphas is not only useful for rejecting 
events produced by the recoiling nuclei from alpha de- 
cays, but is also useful for estimating the amount of con- 
tamination present that may produce other backgrounds. 
In particular, we can estimate the amount of ^-'^^Pb on our 
detectors by detecting alphas produced by the decay of 
one of the daughters of the ^^"Pb decay chain, ^^"Po. In 
the case where the ^-'^"Pb and ^^"Po are in equilibrium, a 
measurement of the alpha rate from the decay of ^^"Po 
yields a beta rate coming from the decay of ^^°Pb and 
•^^°Bi. We have used the measured alpha rates during the 
WIMP search together with simulations of beta contami- 
nants to estimate what fraction of our background is due 
to betas from the ^^°Pb decay chain. 

In the following sections, we describe our simulations 
and measurements of these backgrounds in more detail. 
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FIG. 11: Plot of ionization energy versus recoil energy for 
events during the WIMP search after the initiation of an old- 
air purge between the outer copper can and the surrounding 
mu-metal shield. Alphas (black dots) interacting in the de- 
tectors are easily identified. Nonlinearity in the phonon en- 
ergy response at high energies is evident in the gammas (gray 
dots). This nonlinearity results in the recoil energy of the 
alphas being underestimated by about a factor of two. 
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A. Neutron Background 

Neutrons arise primarily from natural radioactivity, 
such as (a, n) reactions from uranium and thorium nu- 
clei and their decay products, and from muon interac- 
tions with nuclei in the apparatus and surrounding ma- 
terials. While neutrons produced by radioactivity in the 
rock are abundant, they have energies of only a few MeV 
and therefore have a large cross-section on protons. The 
flux of such neutrons is attenuated by a factor of ~ 10^ 
by the 50-cm-thick polyethylene shield surrounding the 
detectors. This attenuation makes the contribution from 
this background negligible for the data taken thus far, 
and even for the entire planned CDMS-II exposure. Cos- 
mogenic neutrons produced by muons entering the shield 
are vetoed with high efficiency on a per-event basis by the 
scintillator panels surrounding the shield (see Sec. Ill C|l . 

The dominant source of any residual neutron-induced 
background that might be observed during the later 
course of CDMS-II is expected to arise from muon inter- 
actions in the surrounding rock cavern. Such neutrons 
can have energies well above 50 MeV, which is the en- 
ergy at which the polyethylene moderator begins to lose 
effectiveness. These fast neutrons, many of which will 
not be detected in the present scintillator muon shield, 
can penetrate the polyethylene and interact, for exam- 
ple, with a Pb nucleus and produce low-energy secondary 
neutrons that could recoil from nuclei in the cryogenic de- 
tectors. In the past !ll'20|, we have measured the rate 
of unvetoed multiple-scatter events due to these "punch- 
through" neutrons to estimate the expected rate due to 
single scatters. In the data we have taken so far at 
Soudan, we did not detect any multiple-scatter events in 
anti-coincidence with the veto, and so we presently rely 
on Monte Carlo simulations to estimate this background. 
In addition, we simulate veto-coincident neutrons which 
are compared with data and serve as a cross-check on the 
anti-coincident population. We are approaching the ex- 
pected required exposure to see such events, but to date 
we have not detected any veto-coincident neutrons. 

We carry out simulations of the fast neutrons produced 
in the rock lisll b v using a pubHshed spectrum 0, ^ 
in a GEANT Monte Carlo program which includes 
the appropriate libraries for hadron cascades. Our pa- 
rameterization of this spectrum, which is dominated by 
hadron cascades caused by primary muon-induced nu- 
clear spallation and by electromagnetic processes, is given 
by: 

dN _ J 6.05 e(~^/'^^) for 50 < E < 200 MeV 

~ |e(-^/25") for i;> 200 MeV ^' 

where E is the neutron energy. 

The normalization is based on the neutron yield in 
rock at Soudan, which has been measured as (3.3 ± 
1.0) X 10"'' n/(£igcm~^) [i^. This yield is consistent 
with FLUKA |43| simulations of the neutron yield of 
3.56 X lO-'' n/(/^gcm-2) at 270 GeV H scaled to 



2.96 X 10~^ n/(/j,gcm~^) at 210 GeV, the mean energy 
of downward going muons at Soudan The scaling 

is based on the mean muon-energy dependence of the 
fast neutron yield, first pointed out by |50l | and since 
confirmed by a variety of measurements and simulations 
(see |5l| and references therein) . The results of this sim- 
ulation, which model the full geometry of the CDMS- 
II apparatus in the configuration of the data set and 
WIMP-search exposure reported in this paper, predict 
0.051 ±0.024 events in the Ge detectors and 0.024±0.011 
events in the Si detectors from neutrons produced in the 
rock. GEANT 4 and FLUKA simulations of muon inter- 
actions in the rock cavern predict that, of these events, 
at least 60% will be coincident with activity in the veto. 

We also carry out simulations of the neutrons pro- 
duced by muons in the lead shield and copper cryo- 
stat |43| . which represent the majority of material in- 
ternal to the scintillator veto. In addition to the high- 
energy component described by Eq. Q , we simulate the 
low-energy component which proceeds through the gi- 
ant dipole resonance, yielding neutrons with the evapo- 
ration spectrum used to describe muon capture processes 
at low depths. The lower-energy component of the pri- 
mary spallation spectrum can thus be described by the 
functional form |52|: 

dN _ j 0.812 E^/^^ g(-B/i.22) f^j. ^ ^ 4 5 j^^y 

'dE ~ jo.OlS e(-^/9) for 4.5 < E < 50 MeV. 

(8) 

This low-energy component is important for neutrons 
produced in the shielding because there is only 10 cm of 
polyethylene internal to the lead and none inside the cop- 
per. The normalization of the neutron yield is inferred 
from measurements in lead for a mean muon energy of 
110 GeV [53I and scaled with mean muon en- 
ergy, as before. The results of these simulations predict 
that neutrons produced in the shield will yield 1.94±0.44 
events in the Ge detectors and 0.89 ± 0.18 in the Si de- 
tectors for the exposure reported in this paper. Since 
the muons producing these neutrons must pass through 
the scintillator, all of these events would be vetoed. We 
find no such nuclear recoils, which is consistent at 6% 
confidence with the predictions of our simulations, indi- 
cating that the neutron background may be less than our 
conservative estimate. 



B. Gamma Backgrounds 

Simulations of the detector response (see Sec. IIV Cp 
show that the vast majority (>99.99%) of gammas in the 
detector are easily discriminated from WIMPs. We can 
identify the different components of our gamma back- 
ground by comparing our WIMP-search gamma spec- 
trum with that predicted by GEANT 13 Monte Carlo 
simulations of our expected gamma background. The 
gamma background in the 5 keV to 3 MeV region for 
39 days of livetime is shown in Fig. ^1 For energies 
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greater than 2.6 MeV, the rates at Soudan are consider- 
ably lower than those of SUF HHH. This is due to the 
absence of high-energy gamma cascades caused by neu- 
tron capture in the materials surrounding and within the 
Ge and Si ZIP detectors. From Zl to Z6 we observe a to- 
tal of 0, 1, 4, 7, 3 and 3 events, respectively, in the energy 
region 2.7-8.0 MeV. All these events are multiple scatters 
depositing energy in the outer charge electrode. In SUF, 
after the muon veto cut, we observed about 240 events 
per detector for a similar exposure of 39 days and the 
same energy interval. These were due to muon-induced 
neutron activation of the detector materials and are sig- 
nificantly reduced at Soudan. Below 2.6 MeV, which is 
the endpoint of natural radioactivity originating from the 
^^^U and ^'^^Th chains, the gamma background is higher 
by about a factor of 2 compared to SUF. This factor is, 
a priori, consistent with the lack of lead shielding inside 
the Soudan icebox which was present for runs at SUF. 




Ionization Energy (keVee) 



FIG. 12: Measured gamma spectra in the coadded Ge detec- 
tors along with Monte Carlo simulations of the most likely 
contributions: U/Th/K on the inner polyethylene (dotted), 
U/Th/K on the Cu cans (thin gray), Rn outside the mu-metal 
shield (thick black), and the sum of all these contributions 
(thick gray). 

We have simulated the contributions to the back- 
ground from the U/Th/K content of the Cu cans, of the 
inner and outer polyethylene shields, and from ^^"Pb in 
the inner lead shield. These materials had been previ- 
ously screened with HPGe spectrometers, and the results 
were used to normalize the Monte Carlo simulations. In 
the Cu used to build the icebox cans, U/Th/K contami- 
nations of less than 0.5 ppb were measured. For the inner 
and outer polyethylene shields, the measured U/Th con- 
tents were 30ppt and 1 ppb, respectively. No ^°K was 
detected, the upper limit being 0.25 ppm. The upper 
limit on ^^"Pb in the inner lead shield was 0.2 ppb. The 



sum of these normalized contributions yields a gamma 
background smaller by about a factor of four than the 
observed one. 

A likely additional source of gammas, especially the 
observed ^^"^Bi line at 609 keV in the Ge detectors, is 
the decay of radon outside the currently purged volume. 
To keep radon gas outside the space directly surround- 
ing the icebox, we purge this volume with air which has 
been stored long enough in metal cylinders to allow the 
radon to decay (T1/2 = 3.2 days) and the daughters to 
plate out. Radon can still penetrate and diffuse into the 
inner polyethylene shield, however, and Monte Carlo sim- 
ulations show that about 35 Bq/ m^^ of radon gas outside 
the purged volume is sufficient to explain most of the re- 
maining gamma background (see Fig. ll2|l . The measured 
radon levels in the Soudan laboratory in the CDMS ex- 
perimental room show a mean of about 500 Bq/m"^ , with 
large seasonal variations (up to 700Bq/m"^ in summer 
and down to 200 Bq/ m'' in winter) . 

C. Beta Backgrounds 

The dominant background for this analysis comes from 
electron recoils within 10 /im of the detector surface. This 
region is referred to as the dead layer. As described in 
Sec. HID II these events have suppressed ionization col- 
lection and could be mistaken for nuclear recoils. Sur- 
face events come from three sources: gammas interacting 
within the first 10 /im of the surface, low-energy electrons 
ejected from nearby material by high-energy x-rays, and 
electrons produced by radioactive beta decays from sur- 
face contamination. The first two sources are related to 
the ambient gamma background and are negligible for 
the background at Soudan. The third source, however, 
can limit the sensitivity of the experiment. 



1. Detector Response to Betas 

We characterize the detector response to contamina- 
tion betas by using the surface events produced in the 
very large gamma calibration runs. To simulate betas 
produced by calibration gammas, we used a simplified 
geometry of CDMS-II Tower 1 where six ZIPs were ver- 
tically stacked and surrounded by a thin copper cylinder. 
Simulated source photons originated isotropically from a 
spherical surface outside of the copper and with energies 
evenly distributed between 0-1 MeV. For each of these 
events, a weight was assigned based upon the initial en- 
ergy of the source particle for that event, such that the 
detector recoil spectrum calculated by the Monte Carlo 
matched that observed in the gamma calibration. We 
simulated ^^'^Ba and ^°Co gamma calibration data. In 
addition to the gamma calibrations, simulations of possi- 
ble beta-emitter ZIP contaminants were conducted with 
sources of ''"K, ^^"Pb, and "'^^C, where each source is dis- 
tributed uniformly between the pair of Ge ZIPs Z2 and 
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Z3. Each event, as described by the Monte Carlo sim- 
ulation, is reported as a list of positions and energy de- 
positions for individual scatters or "hits." An effective 
single position for each event is formed by taking the 
energy-weighted mean position of the hits. 

In order to compare the results from simulations of 
beta emitters with actual data, we need a method for es- 
timating the detector ionization yield of a Monte Carlo 
event. It has been shown 17] that the dead layer ef- 
fect, which leads to incomplete charge collection near 
the surface of the ZIP detectors, can be described by 
a function relating depth to charge collection efficiency 
(see Fig. I13|l . Combining that function with the Monte 
Carlo event depth information allows us to model the 
yield response of the ZIP detectors to simulated surface 
events. This method of simulating betas has been shown 
to accurately reproduce the detector response to a ^°^Cd 
source |5^. 
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FIG. 13: Charge collection efficiency as a function of depth 
from the surface of a Germanium ZIP at a 3.0 V ionization 
bias, as derived from Monte Carlo simulations and analysis of 
^°^Cd data. 

Those leakage events with a signal in more than one 
detector can be excluded as WIMP candidates because 
WIMPs are not expected to scatter in multiple ZIPs. 
However, leakage events with a signal in just one detector 
do remain a problem. 

The results from the depth distribution analysis of 
the ^^■^Ba calibration simulation (see Fig. I14() show that 
single-scatters have a flat depth profile extending into the 
10 fim dead layer. We therefore expect very large Ba 
calibrations to have a low-yield tail coming from single- 
scatters within this region of suppressed charge collec- 
tion. These events can be used to characterize the de- 
tector response to interactions in the dead layer. In the 
region 15 fj,m from the surface, there is a factor of 3 en- 



hancement in nearest-neighbor double events. Beta-Beta 
Double events (BBD), nearest- neighbor double events 
that have low-yield events in both neighboring detectors, 
show high preference for the surface with an exponential 
decrease and 1/e folding scale of ^ 10 /im. 

Simulations of a beta-emitter (see Fig. I14|l show 
that single-scatters from contamination on the detector 
surface have a depth distribution very similar to that of 
^^•^Ba BBD events, suggesting that we could improve our 
characterization of the detector response to beta-emitter 
singles through BBD events from gamma calibrations. 

We can use the function relating the charge collection 
efficiency to depth to estimate what fraction of contami- 
nation betas we expect to have measured ionization that 
is comparable to that of a nuclear recoil. The simulations 
indicate that ~85% of the betas will have yields below 
the electron- recoil band, shown in Fig. of which ~ 25% 
may be misidentified as nuclear recoils. 
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FIG. 14: Depth distributions for simulated surface events cor- 
responding to ^^^Ba single scatters (dashed light gray), ^^^Ba 
nearest neighbor doubles (solid light gray), ^'''^Ba beta-beta 
doubles (dark gray) and ""^K singles (black). The ^^^Ba simu- 
lation results are normalized to all events. The '"'K simulation 
results have been scaled down by a factor of ~100 for easier 
comparison with the results from the ^^^Ba simulations. The 
horizontal thin black line indicates the expected distribution 
for events interacting uniformly throughout the entire detec- 
tor. 

After initiating the old-air purge at Soudan, the beta 
rates in the ZIP detectors decreased by a factor of 1.5- 
2.5. By relating the decrements in both gamma and beta 
rates, we have estimated that surface-electron recoils pro- 
duced by external photons contribute about 40-50% of 
our beta background. However, simulations of ^■^^Ba and 
^''Co sources and comparison with calibration data in- 
dicate that only about 2% of these surface-electron re- 
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coils will have their ionization sufficiently suppressed that 
they would appear as nuclear recoils (prior to any timing 
cuts), implying that this background is negligible for the 
CDMS-II experiment. 



2. Detector Surface Contamination 

The remainder of the beta rate is most probably caused 
by long-lived beta emitters contaminating the detector 
surfaces. A likely candidate is ^^°Pb, which is a long- 
lived (Ti/2 = 22.3 years) radioisotope in the decay chain 
of airborne ^■^^Rn. Our detectors are exposed to radon 
during fabrication, mounting, and testing. We have mea- 
sured the ambient radon levels at all the locations where 
detectors travel and have taken precautions to store them 
in nitrogen-gas purged cabinets or evacuated vessels for 
most of the time. 

The ^-'^"Pb /3~ contamination decays to ^^"Bi, with an 
endpoint of 63keV. ^^°Bi decays to ^^°Po, which then 
decays with a half-life of 138 days to the stable ^°^Pb, 
emitting an a particle with an energy of 5.3 MeV. As 
shown in Fig. ^2 the alpha decays from the decay of 
^^°Po produce a cluster of events in the detectors above 
1 MeV in phonon recoil energy and with a significantly 
lower ionization yield than electron recoils in the bulk. 
Prior to this WIMP-search run the detectors had not 
been exposed to ambient air for about two years, well 
exceeding the 138-day half-life of ^^°Po, and leaving the 
^^°Pb in equilibrium with ^^"Po. Measuring the alpha 
rate thus translates directly into a total beta rate due to 
radon contamination. 

The first column of TableCllists the total number of al- 
phas detected over the course of the entire WIMP search. 
Our analysis rejects events for which there is a signifi- 
cant amount of charge collected by the outer electrode 
(see Sec. lVB4)l . so only alphas that produced signals in 
just the inner electrode are considered. We will assume 
that we have 100% detection efficiency for alphas and 
that all of the alphas come from ^^"Po decays which are 
in equilibrium with ^^°Pb. Consideration of the two rel- 
evant beta-emitters suggests that we will detect roughly 
the same number of low-energy betas in our detector as 
alphas. Simulations indicate that ~10% of these low en- 
ergy betas will be single scatters having ionization yields 
comparable to nuclear recoils. The second column of Ta- 
ble m lists the expected number of single-scatter betas 
that will have yields comparable to nuclear recoils under 
these assumptions. This analysis indicates that ~35% 
of the single-scatters in the WIMP-search nuclear-recoil 
band may come from beta-decay of ^^°Pb and ^^"Bi. 

Another possible source of beta contamination in 
our experiment is ^"^C, produced in the atmosphere by 
cosmic-ray spallation. Natural carbon could be intro- 
duced during processing, and traces of carbon are eas- 
ily measured by various materials-science surface screen- 
ing methods (e.g., Auger depth profiling, Rutherford 
Backscattering, and Particle Induced X-ray Emission). 



Detector 


Number of 
alphas under the 
inner electrode 


Expected Number 

of betas in the 
nuclear-recoil band 


Zl 


26 


2.6 


Z2 


20 


2.0 


Z3 


14 


1.4 


Z4 


17 


1.7 


Z5 


31 


3.1 


Z6 


39 


3.9 



TABLE I: Alpha rates and associated beta-rate predictions 
assuming the ^^'^Rn decay chain is in equilibrium after the 
^^"Pb plate-out. The first column lists the number of alphas 
under the inner electrode found in the WIMP-search data. 
The second column lists the expected number of betas that 
are misidentified, assuming that we have a 100% detection 
efficiency for alphas and that all of the alphas are from ^^"Po. 



Results of these analyses indicate that about 2.5 mono- 
layers of C are present on the ZIP surfaces, consistent 
with exposure to air, but there is no evidence of buried 
carbon. This yields a predicted ^''C beta-event rate of 
0.3 betas/day per ZIP, with a 156 keV endpoint, or less 
than 10% of the total observed 15-45 keV beta rate. 



V. ANALYSIS 

With the WIMP-search data in-hand, along-with cal- 
ibration and simulations of the expected backgrounds, 
we will now consider in detail the pulse-shape analy- 
sis, cut-setting, efficiency estimation, leakages, and sys- 
tematics. Note that during this whole analysis process 
the WIMP-search data itself is blinded, with a gener- 
ous nuclear-recoil band defined such that any possible 
WIMP-candidate events are not present until "the box is 
opened" so that no human bias enters the construction 
of the data cuts. 



A. Event Reconstruction 

1. Waveform Analysis 

For the analysis of the ionization waveforms, we first 
construct detector-specific templates, both for the pri- 
mary pulse and for crosstalk between the inner and outer 
electrodes. We restrict the template construction to us- 
ing pulses in the 10-100 keV energy range in order to min- 
imize non-linear effects in this energy range of interest for 
a candidate WIMP-recoil event. The templates are con- 
structed by averaging a number of ionization pulses from 
in situ calibrations with external ^'^^Ba photon sources. 
Each of these selected ionization pulses is for an event 
that produced ionization that was collected in only one 
of the two electrodes. We use two algorithms to estimate 
the amplitude and time offset of an ionization pulse. In 
the optimal-filter algorithm [E^, construct a fil- 
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ter in the frequency domain from the pulse template and 
the noise spectrum for a given charge channel. This fil- 
ter weights frequency bins based on their signal-to-noise. 
Convolving the filter with each ionization pulse yields 
an estimate of the amplitude and time offset relative to 
trigger for each pulse. In the second algorithm, we use 
a time-domain minimization of the template for each 
pulse to estimate the amplitude. The quantity we mini- 
mize is 



Peak Quadrant 
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ty 
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tc — to 
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tA — to 


tc — to 



TABLE II: Definition of tx and ty parameters where 
ti (i=A,B,C,D) is the time the pulse in the i**^ channel reaches 
20% of its peak. 



N 

E 

1=1 



(9) 



where Vi are the [N — 2048) digitized data samples, 
Si is the pulse-shape template, Vb is the fitted pulse 
amplitude, and a is the rms noise per sample. The 
time-domain algorithm is "non-optimal," as it assumes 
that the frequency dependence of the noise can be dis- 
carded in the amplitude and time offset estimation pro- 
cess. This algorithm is of interest, however, for high- 
energy (1 MeVee) events with pulses that saturate (ex- 
ceed the range of) the digitizers. The time-domain fit 
is applied only to the unsaturated digitizer samples (the 
rising and falling edges of the pulse), making it a bet- 
ter estimator than the optimal-filter algorithm for these 
high-energy events. For both algorithms, we analyze the 
inner and outer electrode signals simultaneously to ac- 
count for crosstalk between the two channels. The time- 
domain fit was inadvertently applied to some of the low- 
energy events in our initial analysis. At lower energies, 
this fit is more sensitive to non-white noise in the charge 
channels, resulting in a slightly worse energy resolution 
than was obtained using the optimal-filter algorithm. 

For the analysis of the phonon waveforms, we use a 
single phonon-pulse template for each detector type (Ge 
or Si). The templates have the functional form of a dou- 
ble exponential with risetime ~30 (15 ^ls for Si) and 
falltime '^300 fis (150 /is for Si). For each phonon chan- 
nel, we determine the area under the phonon pulse, and 
hence the phonon energy, using two different algorithms 
for different energy regimes. At relatively low energies 
(< 100 keV) we use optimal filtering, where, as with the 
ionization signal, a filter generated from the pulse tem- 
plate and the measured noise is used to estimate the pulse 
area. At higher energies, where physical saturation of 
the phonon sensors leads to large variations of the pulse 
shape, integration of the pulse, after subtraction of the 
average pre-pulse baseline and high-pass filtering, pro- 
vides a more accurate estimate for the area of the pulse. 

For each phonon pulse, we also estimate the times at 
which the pulse reaches 10%, 20%, and 40% of the peak 
along the rising edge. From these phonon pulse parame- 
ters, we construct additional quantities which are useful 
for event position reconstruction. The phonon "delay" 
time ti (z =A, B, C, or D for each channel) is the differ- 
ence between the 20% time of the phonon pulse and the 
similarly calculated 20% time of the ionization pulse. We 
define the phonon "risetime" as the difference between 



the 10% and 40% phonon times. For each detector, the 
"peak delay" t and the "peak risetime" t are defined as 
the delay and risetime of the channel that has the most 
energy, which we call the "peak" sensor. As shown in 
Sec. IVB 91 both the peak delay and the peak risetime 
provide excellent rejection of surface events. 



2. Position Reconstruction 

It is possible to combine the phonon parameters to ex- 
tract information regarding the interaction location for 
any given event. We define the x and y axes as shown in 
Fig. 0] The z-axis corresponds to the orthogonal coordi- 
nate, i.e. depth into the crystal. We calculate x-y delay 
coordinates tx , ty from the relative phonon delays in the 
peak sensor and its two neighbors according to Table ITTl 
Figure [Tsl shows that these coordinates are related to the 
physical location of each event. Under the approximation 
that each sensor is a point far from the interaction region, 
these formulas would yield accurate physical x-y coordi- 
nates as defined in Fig. 01 However, the finite size of the 
sensors results in a non-linear mapping between the pa- 
rameters tx , ty and the physical position. As indicated 
in Fig. 1151 the expected azimuthal symmetry is broken 
by the sharp increase in sensitivity near the boundaries 
between the sensors. The "barrel" shape of the event 
distribution suggests that the radial dimension may not 
be monotonic (see Fig. I17() . Modeling of the dynamics of 
the phonon system to obtain a more linear reconstruction 
from the timing parameters 59| has been achieved, but 
for simplicity, we use the parameters described above for 
our current analysis. 

A second nonlinear measure of an event's x-y position 
is determined from the relative partitioning of energy 
among the four phonon channels. These "phonon par- 
tition" position variables are defined by 
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jPc + Pp) - (Pa + Pb) 
EP^ 



and (10) 

{Pa + Pd) - {Pb + Pc) 



Py 



EP^ 



where Pi {i = A, B, C, D) is the phonon energy measured 
in the ith channel. The layout of the QETs produces a 
square (or "box") distribution of the events in the de- 
tector, as shown in Fig. 1161 There are two important 
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FIG. 15: Plot of the time delay parameters ty versus tx for 
detector Z5 with ^^^Ba calibration data. The barrel shape of 
the plot indicates that this event reconstruction is non-linear. 



features of this distribution. First, the corners of the 
box correspond to the centers of each detector quadrant, 
where we expect to find the maximum partitioning of 
the phonon energy. We note that for any event, enough 
phonon energy is distributed to all of the channels such 
that there is an upper bound for Px,y of ~ 0.5. Second, 
the plot folds over (see Fig. [T7|) with events near the outer 
edge of the detector appearing closer to the center of the 
plot. 
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FIG. 16: Plot of the phonon partitioning parameter Py versus 
Px for detector Z5 with ^'^^Ba calibration data. Py and Px 
denote relative partitioning of the phonon energy as defined 
by Eq. itTUl . 



The non-linearity inherent in both measures of event 
position is evident in the sparsity of events in the cen- 
tral regions of the two plots, and the over-crowding to- 
wards the perimeter (where both measures lose sensitiv- 
ity). The "button-hook" feature shown in Fig. 1171 illus- 
trates that both the phonon-timing and phonon-energy- 
partitioning positions are not single- valued. However, 
the combination of both the phonon-timing and phonon- 
energy-partitioning variables does allow accurate param- 
eterization of an event according to its physical x-y posi- 
tion. This method does not reconstruct an event's depth 
(z position), but instead allows us to account for position 
variation of the peak risetime and peak delay parameters 
so that they can be used to reject surface events (see 
Sec. rvTl . 

For every event, this method defines a three- 
dimensional vector 



Q = {Px,Py,R/Ro), 



(11) 



where R 



X 



ty is the delay "radius," and Rq —22.5 



IIS in Ge (12.25 /is in Si). We chose Rq so that the nor- 
malized radius, R/Rq, would have a maximum value that 
is comparable with the maxima of Px.y ■ This vector pro- 
vides a unique mapping for the x-y position of an event 
as shown in Fig. El 
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FIG. 17: Plot of delay "radius,' 
phonon partitioning parameter Px for 
Events shown are for detector Z5 when exposed to a 
source. 



-0.3 < Py < -0.2. 
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3. Position Dependence of the Detector Response 

The phonon pulse shapes depend on the event location 
within the detector. There are two reasons for this vari- 
ation. The first is that the superconducting transition 
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temperature of the tungsten TESs is not perfectly uni- 
form across the face of a detector. Since we bias the TESs 
in parallel, TESs with different transition temperatures 
are biased at slightly different points within their tran- 
sitions, and this varies their response. Detector Zl has 
a very large variation in its TES transition temperatures 
which dominates the position dependence of the phonon 
response. The other detectors in this tower had their 
transition temperatures tuned ^3 to be sufficiently uni- 
form so that this contribution to the position dependence 
of the phonon response is negligible. The second reason 
why the phonon-pulsc shapes vary is that the physical 
arrival time of phonons at the QETs depends on posi- 
tion. For events near the center of a detector, a large 
number of the phonons detected have their last scatter- 
ing point close to the initial interaction location. For 
events closer to the edge of the QET array, more of the 
phonons detected have scattered from the crystal surface. 

The varying phonon pulse shape leads to a position 
variation of our extracted event parameters. Figures ITHl 
1191 and 1201 show that the peak delay, peak risetime, and 
the phonon amplitude parameters vary with position. 
Such variation, if left uncorrected, would degrade the en- 
ergy resolution and event-type discrimination based upon 
the ionization yield. In addition, variation in the tim- 
ing parameters would prevent the use of a single cut to 
discriminate surface events from events in the detector 
bulk. To remedy this, we apply position-dependent cor- 
rections '61^ to make each parameter independent of the 
x-y position of the event. 

First, we correct for a non-linearity of the phonon en- 
ergy response arising from the saturation of the tungsten 
TESs. The correction has the form 



Pt 



A(e 



Pr/A 



1) 



(12) 



where Pt is the total phonon energy and A is determined 
empirically based on the fact that the ionization energy 
approaches linearity at low energies (see Section Fill C|l . 
We correct for a weak energy dependence of the phonon 
timing parameters by using 

(13) 



a + bP^' 

where r and r correspond to the energy-corrected and 
original timing parameters and Pt is the total (unlin- 
earized) phonon energy. We empirically determine the 
constants, a, b, and c such that the corrected timing pa- 
rameter, T, is linear and has the following normalization: 
6 fis (3 lis for Si) for the peak delay, and 12 /is (6 /is 
in Si) for the peak risetime. Finally, we use a lookup- 
table algorithm (described below) to remove the position 
dependence of these energy-corrected quantities. 

A position correction for each event is obtained us- 
ing a lookup table specific to each detector. Each 
lookup table consists of approximately 12,000 single- 
scatter inner-electrode electron recoils with energies be- 
tween 10-200 keV. These electron recoils come from cal- 
ibrations with external photon sources. For each event. 



i, in the lookup table, we determine a "distance," pij, to 
any other event, j, in the table according to the following 
metric: 



\ 



3 



k=l 



(14) 



where Q is as defined by Eq. (|11() . For each event in the 
lookup table, we determine the N nearest events in the 
table. The value of N is the same for every event in a 
given detector. N must be sufficiently small to remove 
the variation due to the x-y position dependence while 
being sufficiently large to prevent statistical variations 
in the correction from contributing significantly to the 
position-independent resolution. For these lookup tables, 
we chose N to be around 80. 

The correction value for the phonon energy of a par- 
ticular event in the lookup table is the average of the 
ionization yield for the event's N nearest neighbors. The 
correction values for the timing parameters are the aver- 
age of that timing parameter's values for the event's N 
nearest neighbors divided by that parameter's normal- 
ization. Events that are not part of the lookup table use 
the correction values for the nearest event in the lookup 
table. 
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FIG. 18: Plot of Peak Delay versus phonon partitioning pa- 
rameter Px for the slice —0.3 < Py < —0.2. Gray dots are 
uncorrected distributions and black dots are corrected distri- 
butions. The corrected Peak Delay has been rescaled to 3 fis 
for easier comparison. 



B. Cut Definitions 

In this Section we describe all the cuts of our analysis 
and estimate their efficiencies. Nearly all efficiencies are 
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FIG. 19: Plot of Peak Risetime versus phonon partitioning 
parameter Px for the slice —0.3 < Py < —0.2. Gray dots are 
uncorrected distributions and black dots are corrected distri- 
butions. The corrected Peak Risetime has been rescaled to 
9.5 fis for easier comparison. 



we inadvertently used the time-based fit to determine 
the amplitude of some of the ionization pulses. To ac- 
count for this, we calculated two sets of efficiencies: one 
for the initial implementation of the cuts, the other for 
the current analysis which applies the cuts as intended. 
Only the fiducial-volume cut for detector Z5 was found 
to have a significantly different efficiency. We will show 
the reduced efficiency in Figs. I^Hl and 123 and describe 
both estimates in Section FV B 41 Tables IVl and IVll sum- 
marize the effect of the ten individual data analysis cuts 
described in detail below. 



1. Data- Quality Cut 

A data-quality cut removes all non-optimal data sets 
(or parts of data-sets). This cut excludes data-sets with 
known problems (such as non-operational channels), data 
sets that failed in the off-line diagnostics (usually due to 
increased noise), and events that were triggered by noise 
bursts. Figure [3 shows that, over the course of the entire 
run, 52.6 live-days of WIMP-search data pass the data- 
quality cut. Less than 5% of the data fail the data-quality 
cut. 
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2. Phonon Pre-trigger Cut 

The phonon pre-trigger cut rejects events for which the 
pre-trigger part of a phonon trace (about 400 /its long) is 
unusually "noisy," more specifically, for which the stan- 
dard deviation of the ADC values for the first 400 0.8 /is 
time bins is large. This cut complements the ionization 
cut described below for purposes of eliminating pile- 
up events and noise-triggered events. For the duration of 
this run, the phonon pre-trigger cut was set at 5tT, where 
(T is the standard deviation of a Gaussian fit performed 
separately for each phonon channel of each detector using 
data from ^'^^Ba gamma calibrations. 

The efficiency of the phonon pre-trigger cut is calcu- 
lated as the fraction of WIMP-search events passing the 
cut. The efficiency is better than 99.999%. 
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FIG. 20: Plot of Ionization Yield versus Px for the slice 
—0.3 < < —0.2. Gray dots are uncorrected distribu- 
tions and black dots are corrected distributions. The cor- 
rected yield has been rescaled to 0.95 for easier comparison. 



estimated as a function of the event's recoil energy in the 
following energy bins: 5-8, 8-12, 12-16, 16-20, 20-30, 
30-40, 40-50, 50-60, 60-70, and 70-100 keV. As men- 
tioned earlier, all cuts were finalized prior to inspecting 
the WIMP-search data set. For the initial analysis |23l |. 



3. Ionization Cut 

Events with anomalous ionization pulse shapes are re- 
jected by a cut on of the optimum-filter fit. Such 
events can be from pile-up (in which more than one 
pulse takes place in the same 1.6 ms trace) or from noise 
glitches. Since the start times of the optimal filter's tem- 
plates can vary only by integer numbers of time bins, 
the fits are slightly inaccurate for events that start be- 
tween two bins. For high energy pulses, this inaccuracy 
is large compared to the noise, so this sub-digitizer-bin 
time-jitter causes some pulses to have worse fits. To keep 
these events, the cut is implemented as a quadratic func- 
tion of the ionization energy, as illustrated in Fig. [2] for 



22 



detector Z3. The cut is defined using the ^'^^Ba cahbra- 
tion data. 

Figure 12] shows the existence of two bands. As indi- 
cated by Fig. 122 the upper band corresponds to events 
during times of increased noise produced by a heater on 
the cryostat. This heater intermittently introduces ad- 
ditional noise at a few hundred kHz in the ionization 
channels and causes worse fits for the events occurring 
while it is on (approximately 5% of the total). Though 
this additional noise did not impact the energy resolution 
of the optimal-filter algorithm, we decided to reject these 
events in order to simplify the cut. 
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FIG. 21: Reduced (for 4092 de grees of freedom) of the ion- 
ization pulse as given by the optimal-filter algorithm plotted 
against ionization energy (electron equivalent) from a ^^^Ba 
calibration for detector Z3. Only events beneath the solid line 
pass the cut. The faint upper band consists of events taken 
during the ~ 5% of the time that a noisy heater was on; see 

Fig.m 
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FIG. 22: Plot of the reduced (for 4092 degrees of freedom) 
of the ionization pulse from the optimal-filter algorithm versus 
time of day for a ^^^Ba calibration of detector Z3. Events are 
between 10-60 keVee. The events in the upper band of Fig. 1211 
are clearly localized in time. These excursions are correlated 
with the time a noisy heater was active. 
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The efhciency of the ionization cut is estimated 
using the WIMP-search electron-recoil events. We use 
WIMP-search data rather than the calibration data be- 
cause we have observed small variations in the ionization 
parameter among the WIMP-search data sets. The 
efficiency is estimated as a function of ionization energy. 
We calculate the efficiency as a function of recoil energy 
for nuclear-recoil events using the observed relationship 
between ionization and recoil energy for neutrons from 
calibration with external sources. As an example, Fig. 1231 
shows the efficiency of these two cuts on the nuclear- recoil 
events for the detector Z3. 
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FIG. 23: Efficiency of the ionization cut for nuclear-recoil 
candidate events (see Sec. lV B'St in the Ge detector Z3. Note 
that the plot is versus recoil energy and not ionization. 



4- Fiducial- Volume Cut 

As described in Sec lIIPTl we use the amplitude of the 
pulse in the outer ionization electrode, or "guard ring," 
to reject events that take place close to the edge of the 
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crystal. The cut is defined in the Qo-Qi plane, where 
Qo and Qi stand for ionization energy in the outer and 
inner ionization channels, respectively. ^'^^Ba calibration 
data are used to fit the Qo distribution to a Gaussian in 
several bins of Qi energy. The events passing the cut lie 
in a band centered around Qo ~ 0, defined by 



Time-based Fit 



fj. — 3a < Qo < H 



./25, 



(15) 



where fi and a are the (Qi-dependent) mean and stan- 
dard deviation determined by the Gaussian fits. The last 
term is introduced in order to reduce the impact of the 
outer-electrode noise on the high-energy inner-electrode 
events. Figure illustrates this cut for the Ge detector 
Z5. 

Optimal Filter 
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FIG. 24: Ionization collected in the Outer electrode versus 
ionization collected in the Inner electrode using the optimal- 
filter (see Sec. IV A l|l algorithm for Ge detector Z5. Cali- 
bration data taken with an external ^'^'^Ba source are shown. 
Events between the two lines pass the fiducial volume cut. 

The efficiency of this cut for WIMPs should be nearly 
82%, which is the volume fraction covered by the inner 
electrode. The efficiency should be slightly less than 82% 
(especially at high energies) because events near the gap 
between the inner and outer electrodes may have some 
ionization collected in the outer electrode. Additionally, 
noise in the outer electrode will cause some events in 
the inner electrode to fail the fiducial volume cut. To be 
conservative, we used the expected efficiency for neutrons 
(which could be checked more directly) as the efficiency 
for WIMPs as well. The efficiency of the fiducial-volume 
cut for neutrons should be lower than the efficiency for 
WIMPs due to two effects, both of which are related 
to WIMPs' small interaction probability. First, some 
neutrons multiple-scatter within a detector and thereby 
deposit energy in both the inner and outer electrodes. 



50 



> 40 



30 



20 



10 



w 

I 

o 



-10 




-50 50 100 150 200 250 
Inner Electrode Ionization Energy (keVee) 



300 



FIG. 25: Ionization collected in the Outer electrode versus 
ionization collected in the Inner electrode using the time- 
based fitting (see Sec. IV A H algorithm for Ge detector Z5. 
Calibration data taken with an external ^^^Ba source are 
shown. Events between the two lines pass the fiducial vol- 
ume cut. 



Second, due to the shielding of the inner region of the 
detectors by their outer regions and by other detectors 
in the stack, neutrons that single scatter are more likely 
to deposit energy in the outer electrode than would be 
calculated by the volume fraction alone. 

We estimate the efficiency of the outer-ionization- 
electrode cut using a Monte Carlo simulation of the ^^^Cf 
neutron calibration. The simulation determines the re- 
coil energy and ionization energy of each detector hit and 
sums the inner and outer ionization energy of all hits in a 
detector based on the simulated positions of the hits. For 
this simulation, we do not account for any possible shar- 
ing of ionization between the two electrodes. Any hit had 
its full ionization collected by the electrode it was closest 
to. Then we convolve the ionization energy with the ob- 
served noise performance of the ionization channels. We 
characterize the noise of each channel by a Gaussian fit of 
the distribution of ionization-energies for random-trigger 
events. We then estimate directly the fraction of the neu- 
trons that fail the outer-ionization-electrode cut. To es- 
timate the efficiency for the time-based peak-height anal- 
ysis, we simply use the noise performance given by the 
time-based algorithm. Only Z5 had significantly worse 
noise with this algorithm (see Fig. I25() . The increase in 
noise results in the initial implementation of this cut be- 
ing more severe than intended, reducing its efficiency as 
shown in Fig. I^Hl 

Figure [2ni also shows a check of our calculated fiducial- 
volume cut efficiency by comparison with a measurement 
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FIG. 26: Efficiency of the fiducial- volume cut for detector Z5 
(Ge). The black line corresponds to the Monte-Carlo sim- 
ulation of the ■^^^Cf neutron calibration estimate when the 
intended optimal-filter algorithm is used to estimate the ion- 
ization pulse amplitude. The gray line corresponds to the 
estimate when using the time-based fit. The black dots with 
error bars are estimates using the actual nuclear-recoil events 
due to neutrons from the ^^^Cf calibration runs using the 
cuts for the current analysis. The roll-off above 60 keV is 
an artifact due to surface-electron contamination of the data 
set. Only the fiducial- volume cut for this detector (Z5) varied 
substantially (~25% reduction) with the the time-based fit. 



using ^^^Cf neutron calibration events. These events were 
preselected to pass the ionization x^, phonon prc-trigger, 
and ionization threshold cuts and to have ionization 
yields within the 2a nuclear-recoil band (see Sec. lV B'Sll . 
Our efficiency estimate is the fraction of the preselected 
events that pass the outer-ionization-electrode cut. The 
efficiency estimated in this way suffers from electron- 
recoil misidentification in the outer electrode. Electron- 
recoil events that take place close to the edge of the crys- 
tal may not have full charge collection and can leak into 
the nuclear-recoil band, contaminating the sample of true 
neutron events. The efficiency is consequently underes- 
timated, especially at high recoil energies (> 70 keV) 
where the number of neutrons in the neutron calibration 
is relatively low. 



5. Electron-recoil and Nuclear-recoil Bands 

To calculate the electron-recoil bands and the nuclear- 
recoil bands, we first fit Gaussians to the distributions 
of ionization yield (= ionization / recoil energy) for both 
electron-recoil and nuclear-recoil events (from the ^'^^Ba 
and ^^^Cf calibrations, respectively) in several recoil- 
energy bins. The estimated means and standard devia- 



tions are then fitted versus recoil energy. Unless specified 
otherwise, the bands are always taken to be ±2cr. 

We determine the efficiency of the nuclear-recoil-band 
cut using the nuclear-recoil events in the neutron cali- 
bration. We preselect events by imposing the charge x^, 
phonon pre-trigger, charge threshold, fiducial-volume, 
and phonon-timing cuts. Additionally, we only consider 
events in the 4cr nuclear-recoil band, which are expected 
to include almost all of the true nuclear-recoil events. 
The efficiency estimate is the fraction of the preselected 
events that fall in the 2(t nuclear-recoil band. Figure 1771 
shows the efficiency of this cut for Z3. Neutrons which 
multiply scatter in the same detector tend to have lower 
yield than do single scatters with the same total recoil 
energy. This effect should cause our measured nuclear- 
recoil-band cut efficiency to be slightly lower than the 
efficiency for WIMPs. We make no correction at this 
time for this small effect. 
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FIG. 27: Efficiency of the nuclear-recoil-band cut for the Ge 
detector Z3 (see Sec. IV HSU . 



6. Ionization Threshold Cut 

The goal of the ionization threshold cut is to select 
events that have measurable pulses in the charge chan- 
nels. This cut is particularly important at low energies, 
for which an event with only noise in the inner ioniza- 
tion channel could otherwise be mistaken for a low-energy 
nuclear-recoil event. 

This cut is defined using the distribution of ionization 
energy for noise events (obtained by random triggering) 
in the ^^'^Ba calibration data. In particular, we exam- 
ine the distribution of the optimal-filter amplitudes of 
pulse-less traces in the inner charge electrode. We fit a 
Gaussian to this distribution and set the cut at /i -I- 5a, 



where /i and a are the mean and the standard deviation 
estimated in the Gaussian fit, respectively. 

The efficiency of the ionization threshold cut is mea- 
sured with inner-electrode, ^^^Cf calibration events which 
trigger on the phonon signal and fall in the nuclear recoil 
band. The efficiency is the fraction of such events which 
pass the ionization threshold cut. The efhciency for Z3 
is shown in Fig.|2Sl Figure also checks this result by 
comparison with a calculation of the efficiency of a nearly 
equivalent cut in the yield vs recoil-energy plane, assum- 
ing the yield distribution is Gaussian. For events with 
zero outer electrode signal, the ionization cut is equiv- 
alent to a hyperbolic cut in the yield vs recoil-energy 
plane. This alternative efficiency estimation suffers from 
a systematic error because the yield actually depends on 
the sum of the inner and outer charge channels. But 
its results are nonetheless very similar to those obtained 
using real events. 
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FIG. 28: Efficiency of the ionization-tlireshold cut (see 
Sec. IV B 6ll on nuclear-recoil events. Black points with er- 
ror bars indicate the efficiency as determined from data using 
nuclear recoils from the ^^^Cf calibration runs. Diamonds cor- 
respond to the calculated efficiencies assuming that the yield 
is a Gaussian. 



7. Muon-Veto Cut 

We define the muon-veto cut using the time between 
the global trigger and the last hit in the muon veto. 
Events for which this time is long are unlikely to be 
caused by a muon interaction inside the veto. As de- 
scribed in Section III Gl events occurring within 50 /iS 
after veto activity are rejected, resulting in removal of 
> 99.4% of events caused by muons entering the shield- 
ing. The size of this window is set primarily to en- 
sure that low-energy muon-induced events which have 
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Time since last veto hit (|lsec) 



FIG. 29: Histogram of time from last veto hit for Tower 1 
at the Stanford Underground Facility. An excess above the 
5.7-kHz rate expected for uncorrelated events is seen for times 
near zero. The cut at —50 /xs is set in order to remove the ex- 
cess for single scatters with recoil energies between 5-10 keV 
in Ge detectors Z2, Z3, and Z5 (solid dark gray curve); in Zl 
(solid light gray line); and in the two Si detectors Z4 and Z6 
(dashed gray line). 

delayed global triggers (see Sec. Ill E|l are removed. Fig- 
ure[55|shows a histogram of the delays for these detectors 
when operated at the shallow Stanford Underground Fa- 
cility |23| , indicating that setting the cut at 50 \is should 
be sufficient even for energies as low as 5 keV. 

The efficiency of the mtion-veto cut for WIMP events is 
determined by considering the effective deadtime induced 
by each veto hit. Given the 600 Hz rate of the muon veto 
(dominated by ambient photons) and the 50 /iS rejected 
window after each hit, about 3% of the data- acquisition 
time is lost. The efficiency of the muon-veto cut is there- 
fore 97%, independent of the energy deposited in a de- 
tector. 



8. Singles and Multiples Cuts 

To define single-scatter events (events in which only 
one detector was hit), we use the distribution of phonon 
energy for noise events obtained by random triggering. 
In particular, single-scatter events are events in which 
only one detector had a phonon signal larger than 6a of 
this distribution. 

The mtiltiplc scatter events have to satisfy the follow- 
ing criteria: 

• ALL detectors have to pass the data-quality cut, 
the ionization cut, and the phonon pre-trigger 
cut. 
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• At least two detectors have to pass their fiducial- 
volume cut and have recoil energies above 5 keV. 

Note that these two definitions do not account for all 
events. For example, an event with 6 keV and 4 keV recoil 
energies in two detectors would not qualify as a multiple 
scatter (< 5 keV in the second detector), but it may also 
not qualify as a single scatter since both detectors may 
have phonon signals larger than Gcr of the phonon noise. 

We can place a lower bound on the efficiency of our 
singles cut on WIMPs by considering the event rate and 
the duration of the post-trigger time window. The av- 
erage event rate for all the detectors is 0.1 Hz, and the 
post-trigger digitization window is 1.2 ms, giving a lower 
bound of 99.99% for the singles-cut efficiency. 
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9. Phonon Timing Cut 

As discussed earlier in Sec. Ill D 21 the sensitivity of a 
ZIP detector to athermal phonons produced by an inter- 
action provides for background rejection using phonon 
pulse-shape information. In particular, electron recoils 
near the detector surface result in a larger ballistic frac- 
tion than nuclear recoils in the bulk, both from the larger 
ionization yield and the more rapid down-conversion of 
phonons generated in a near-surface event. Figure[Slillus- 
trates these effects in calibration data. Using the phonon 
pulse shape to reject surface-electron recoils is especially 
significant, since these events can have reduced ioniza- 
tion collection, leading to misidentification if using the 
ionization yield alone. 

A detailed beta calibration of a Ge ZIP detector per- 
formed at the UC Berkeley test facility [s^ indicated that 
using two phonon timing parameters, the peak delay and 
the peak risetime (see Sec. IV A"!)) , provides good rejec- 
tion of surface-electron recoils while retaining reasonable 
acceptance of nuclear recoils. 

To define the phonon timing cut, we first determine 
cuts (see Fig. I3(J|I in each of the two parameters for the 
following recoil energy bins: 5-10 keV, 10-20 keV, 20- 
40 keV and 40-100 keV. In each bin, we determine a pair 
of cut values that would reject every event from a subset 
of ^"^^Ba calibrations that was in the 4(7 nuclear-recoil 
band. Using only a portion of the ^'^^Ba calibrations 
allows us to cross check our misidentification estimates 
with an independent set of events (see Sec. IVD|I . For 
each parameter, we then perform a piecewise linear fit to 
the cut values to define a continuous, energy-dependent 
cut (see Fig. I31II . By defining the cut in this way, we 
expect to misidentify a fraction of an event in each of the 
energy bins (see Sec. lVD|) . 

We estimate the efficiency of the phonon-timing cut on 
WIMPs by using neutrons from the ^^^Cf calibrations. 
We preselect events in the 2a nuclear-recoil band passing 
data-quality cuts (charge x'^, phonon pre-trigger, charge 
threshold) and the outer charge-electrode cut. The effi- 
ciency estimate is the fraction of these events that also 
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FIG. 30: Plot of phonon peak delay versus phonon peak rise- 
time (defined in Sec. IV A 11 for Ge detector Z5 in the 20- 
40 keV recoil bin. ^^^Cf neutrons in the 2a nuclear-recoil band 
(gray) and ^^^Ba surface-electron recoils in the Aa nuclear- 
recoil band (black) are shown. The first step in determining 
the timing cuts is to set energy-independent cuts within the 
bin for each of these two parameters so that all of the surface- 
electron recoils are rejected. Events below or to the left of the 
lines illustrating hypothetical cut values would be rejected. 

pass the phonon-timing cut (see Fig.l^. 



10. Analysis Thresholds 

Effective discrimination of electromagnetic back- 
grounds is achieved for energies as low as 5 keV in most 
detectors. However, to be conservative, we set the anal- 
ysis thresholds for most detectors at 10 keV. The excep- 
tions are Zl and Z4, whose surface-event rejection was 
inadequate for recoil energies < 20 keV. For these two 
detectors, we set the analysis threshold at 20 keV. All 
detectors have an upper analysis bound on the recoil en- 
ergy of 100 keV as this is an appropriate upper limit on 
the WIMP-search candidate-event recoil energy of inter- 
est. 



C. Overall Cut EfHciency 

The cut efficiencies described above are multiplied to- 
gether to obtain the overall cut efficiency. Note that the 
efficiency is forced to zero below 10 keV for all detectors, 
and below 20 keV for Zl and Z4, to reflect their analysis 
thresholds. Furthermore, Z6 is not used in this analysis 
because of the ^'^C contamination mentioned in Sec. IIVI 
and neutralization problems mentioned in Sec. Ill DH We 
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FIG. 31: Plots of peak risetime (top) and peak delay (bot- 
tom) versus recoil energy for detector Z5 (Ge). Gray dots 
correspond to ^^^Cf neutrons in the 2a nuclear-recoil band. 
Black dots correspond to ^^^Ba surface-electron recoils in the 
4(T nuclear-recoil band. The piecewise linear functions corre- 
sponding to the cut in each parameter are also shown. Events 
below these lines are cut. For this detector, above ~20 keV 
all of the discrimination comes from the peak delay. 



average over the Ge detectors Zl, Z2, Z3 and Z5 and per- 
form a fit to the efficiency estimate versus recoil energy: a 
straight line above 20 keV and a parabola below 20 keV, 
vifith the slopes matched at 20 keV. Figure shows the 
effect of imposing each individual cut as well as the over- 
all efficiency estimate for this analysis. The discontinuity 
at 20 keV comes from the higher analysis threshold of Zl. 
Integrated over the 10-100 keV bin, the overall efficiency 
yields a net Ge exposure of 19.4 kg-days for a 60 GeV c~^ 
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FIG. 32: Efficiency of the phonon-timing cut for the Ge de- 
tector Z3. 



WIMP. 

For the initial analysis, where we inadvertently used 
the time-based algorithm to determine the ionization en- 
ergy for some of the data, we had to calculate the effi- 
ciencies for both of the algorithms and make a weighted 
average of the two, where the weights were determined 
from the fractions of the WIMP-search data that were 
analyzed by each algorithm (see Table lllT|l . As shown in 
Fig. 1331 the efficiency for this analysis is slightly lower 
than the efficiency for the analysis using just the optimal 
filter algorithm. As discussed in Sec. IVB 41 the reason 
for this loss in efficiency is that the time-based fit leads to 
worse resolution for the outer electrode in Z5. This reso- 
lution loss made the guard ring cut for this detector more 
severe, reducing the efficiency of the cut. The efficiency 
for the Si detector (Z4) was unaffected (see Fig. I34() . 



Detector 


Fraction analyzed by time-based fit 


Zl 


66.1% 


Z2 


57.2% 


Z3 


50.1% 


Z4 


41.0% 


Z5 


51.2% 


Z6 


37.0% 



TABLE III: Fractions (in live-time) of the WIMP-search data 
analyzed by the time-based fitting algorithm, for all six de- 
tectors. The remainder of the data was analyzed by the 
frequency-based optimal-filter algorithm which has better res- 
olution at low energies. 



28 




20 40 60 80 

Recoil Energy (keV) 



100 




20 40 60 

Recoil Energy (keV) 



100 



FIG. 33: Overall WIMP efficiency when cuts are applied. 
Light gray dashed line corresponds to application of data- 
quality, pile-up, and muon-veto cuts and, above 20 keV, over- 
laps with the light gray solid line which corresponds to impos- 
ing the ionization analysis threshold. Dark gray dashed line 
corresponds to application of the fiducial-volume cut. Dark 
gray solid line corresponds to the 2cr nuclear-recoil-band cut. 
The black solid line corresponds to application of the phonon 
timing cuts and is the overall efficiency for our current anal- 
ysis. The black dashed line is the overall efficiency for the 
initial analysis. 



D. Background Rates and Leakages 

The backgrounds for the WIMP-search results pre- 
sented in this paper consist of two classes of interactions 
that could be misidentified as WIMPs: nuclear recoils 
from ordinary particles and electron recoils in the dead 
layer. As discussed in Sec. IIV Al Monte Carlo estimates 
of the punch-through neutron background at Soudan pre- 
dict 0.051 -I-/- 0.024 events in the coadded Ge detectors 
and 0.024 -f/- 0.011 events in the Si detector over the 
course of the entire run. Of these events, more than 
60% will have activity in the veto and will be rejected 
by the muon-veto cut. This gives an expected neutron 
background of 0.018 ± 0.008 in Ge and 0.004 ± 0.002 in 
Z4. Nuclei recoiling from alpha decays are also expected 
to produce negligible background because the efficiency 
for detecting an emitted alpha particle is high. As dis- 
cussed in Sec. IIVI gammas and betas constitute a back- 
ground through interactions in the dead layer where the 
ionization collection can be incomplete. In this section, 
we describe our estimated misidentification of dead-layer 
events. Table IIVI summarizes our expected background 
misidentification. 

We defined the timing cuts to reject all of the events 
in the 4cr nuclear-recoil band coming from a subset of the 



FIG. 34: Overall WIMP efficiency for the Si detector Z4. 
Dashed line corresponds to the efficiency for the initial analy- 
sis and the solid line is the efficiency for the current analysis. 
The two estimates differ by less than 1% since, for this detec- 
tor, there was little difference between using the time-based 
pulse height algorithm and the optimal filter algorithm. 



^^•^Ba calibrations (see Sec. lVB 9)l . Since the cut is set to 
reject all of the events, we expect, on average, to misiden- 
tify one event for every N]^ + 1 events, where 7V2t is the 
number of events in the 4(t nuclear-recoil band that were 
used to determine the cut. It is important that we sepa- 
rately estimate our leakage for each of the energy bins de- 
fined in Sec lVBOl If the surface events produced during 
^^■^Ba calibrations yield distributions of the phonon tim- 
ing parameters that are reasonable approximations to the 
distributions of these timing parameters for the sources 
of beta contamination (see Sec. IIV C|l . we can estimate 
the number of misidentified surface-electron recoils by 
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(16) 



where n2a is the number of events in the WIMP-search 
data passing all but the timing cuts. This leakage esti- 
mate gives an expected surface-electron-recoil misidenti- 
fication of 0.7 ± 0.3 events for Ge and 0.15 ± 0.12 events 
in Si (Z4 only) for our entire WIMP-search data. 

As a crosscheck of our algorithm for estimating 
misidentification, we use Eq. (|16|l to calculate the ex- 
pected leakage for a set of ^^■^Ba calibrations that were 
not used to define cuts. For these calibrations, there are 
a total of 141 events in Ge and 51 in Si (Z4 only) that 
are in the nuclear-recoil band. Substituting these num- 
bers for n2cr (in the appropriate energy bins) in Eq. H16() 
we expect 4.55 ± 1.44 events in Ge and 0.42 ± 0.31 events 
in Si to pass the timing cuts. In the calibration data, 
three events pass the timing cuts in Ge and one passes 
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the timing cuts in Si, consistent with our calculated leak- 
age. This agreement gives us confidence in our method of 
estimating our surface-electron-recoil misidentification. 

We note that there could be an important system- 
atic associated with this calculation. Misidentified events 
from the ^■^^Ba calibrations consist entirely of surface 
events produced by ambient gammas. As discussed in 
Sec. IIVI dead-layer events produced during the WIMP 
search are dominated by beta decay from radioactive 
contamination on the detector surface. Differences in 
the spatial distribution between dead-layer events com- 
ing from contamination and events produced by ambient 
gammas can lead to a systematic error in our leakage 
estimate. We utilize the contamination on the silicon de- 
tector Z6 to estimate the impact of this systematic for 
both Ge and Si detectors. 

Over the course of the entire WIMP search, there were 
150 single-scatter events in the Z6 nuclear- recoil band 
that were anticoincident with the veto. Of those events, 
we expect fewer than 15 to be produced by the ambient 
gammas, indicating that the majority of the events came 
from contamination on the detector itself. The estimate 
of our surface-event misidentification predicts that, out of 
these 150 events, 3.06±1.23 events should pass the timing 
cuts. There were five events that pass the timing cuts, 
which is consistent with our expectation. This agreement 
suggests that the differences between surface events pro- 
duced by contamination and surface events coming from 
ambient gammas do not lead to significant differences in 
discrimination by the timing cuts. 





Ge 


Si 


surface-electron recoils 


0.7 ±0.3 


0.15 ±0.12 


Neutrons 


0.018 ± 0.008 


0.004 ± 0.002 



TABLE IV: Total number of background events appearing in 
the nuclear-recoil band after cuts expected for the 52.6 live- 
days of the Tower 1 WIMP-search run. 



E. Systematics 

Small variations in the detector performance through- 
out the course of the WIMP search lead to a systematic 
uncertainty in our calculation of the cut efficiencies, since 
these efficiency estimates are determined from the brief 
neutron calibrations. We estimate this uncertainty by 
comparing the cut efficiencies on the electron recoils over 
the course of the WIMP-search data with the efficiencies 
on electron recoils during the neutron calibrations. The 
systematic uncertainty can be divided into two parts: one 
component comes from variation in the calibration of the 
ionization and phonon measurements, and one compo- 
nent comes from variation in the detector phonon timing 
response. Despite the use of the lookup-table correction 
described in Sec. lV A'3l there may still be a residual po- 
sition dependence which would correlate the two effects. 



In order to properly account for any such correlations, 
we compare the efficiencies of the two cuts combined. 

Estimating the uncertainty in the timing-cut efficien- 
cies is especially difficult since, as shown in Fig. |S1 the 
timing distribution of the electron recoils is intrinsically 
different from that of the nuclear recoils. To compute 
the uncertainty due to variation in the phonon timing 
response, we first parameterize the electron-recoil tim- 
ing distributions as Gaussians with means and standard 
deviations given by: 

MmMP = ^^cf + AAi, (17) 

and 

<^WIMP = f^cf X (18) 

where /i^^ (wimp) ^'^^ '^cf (wimp) ^^'^ mean and 
standard deviation of the distribution in the neutron cal- 
ibrations (WIMP-search data). Variations between the 
electron recoils of the calibration and those of the WIMP 
search are described by a shift in the means of the distri- 
butions, A/i, and a scaling of the standard deviations, A. 
We then adjust the timing distributions of the neutron 
calibration according to the equation: 

t' = Ax(t-A^^«) + M^; + A^, (19) 

where iJ-cp and cr^^ are the mean and standard devia- 
tion of the nuclear-recoil timing distributions in the ^'"^^Cf 
data. Equation (|19|l simply shifts the mean of the timing 
distributions by A/i and scales the standard deviation by 
A. We estimate the variation of our timing-cut efficiency 
to be the fraction of events that pass the timing cut using 
these adjusted timing parameters. 

The results of this computation of our systematic un- 
certainty indicate that our efficiency estimates are good 
to within 10% in the energy range of 10-100 keV. These 
estimates are consistent with other checks of experi- 
mental stability including comparison of the ionization 
peaks in three separate ^'^^Ba calibrations, observation 
of the 10 keV line from Ga activation throughout the 
WIMP search, analysis of the timing distribution of sur- 
face events over the course of the WIMP search, mon- 
itoring of the bandwith of the TES-readout electronics, 
and comparison of the nuclear-recoil timing distributions 
for the three separate neutron calibrations. 



VI. RESULTS AND DISCUSSION 

In the 52.6 live days for this first WIMP search at 
Soudan, we find one event with a recoil energy of 64 keV 
passing all of the cuts in our current analysis (see Fig.l36|l. 
We note that the timing cuts are much more severe in re- 
jecting gammas compared to nuclear recoils (see Fig.lSl, 
leading to a significant reduction of events in the electron- 
recoil band. Figure ISSl illustrates that our analysis still 
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retains sensitivity to nuclear recoils while rejecting a sig- 
nificant fraction of the electron recoils. For the initial 
analysis where we used the time-based algorithm to de- 
termine the ionization peak height, we found no events 
passing all of our cuts (sec Fig. 13 7() . These two results 
are consistent with each other, given the differences in 
exposure and statistical uncertainty. The observation of 
one event at 64 keV is consistent with our expected (sur- 
face) electron recoil misidentification. Tables and IVII 
show the number of events passing the cuts for the Ge 
and Si detectors, respectively. We note that several new 
advanced analysis techniques (e.g. [l^), which are be- 
yond the scope of this paper, are substantially improving 
the rejection of surface electrons, so that we do not ex- 
pect to be limited by this background for the remaining 
CDMS-II runs. 





Initial 


Current 


All events 


968,680 


968,680 


Not random trigger 


940,619 


940,619 


Analysis thresholds { Sec. IV B 101 


79,655 


79,460 


Singles fSec.lVBSfl 


20,715 


20,907 


Data quality (Sec. IV B 111 


18,852 


19,027 


Pile UD fSects.lVB2landlVB3l 


17,622 


17,793 


Muon veto f Sec.lVB7ll 


17,171 


17,339 


Ionization threshold (Sec. IV B 61 


14,697 


14,835 


Fiducial volume ( Sec. IV B 41 


7,187 


7,615 


Nuclear-recoil band fSec. IV B 5ll 


29 


23 


Phonon timing (Sec.lVBQt 





1 



TABLE V: Breakdown of the events in the Ge WIMP-search 
data as we apply each cut. One event with a recoil of 64 keV 
passes all of the cuts in the current analysis. 
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All events 


968,680 


968,680 
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Sineles fSec.lVB8l 
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3,163 
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Fiducial volume (Sec. IV B41 
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Nuclear-recoil band (Sec. IV B 5ll 


14 
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Phonon timing (Sec. IVB91 









TABLE VI: Breakdown of the events in the Si WIMP-search 
data as we apply each cut. 

Under the assumptions of a standard galactic halo (as 
described in these data yield the lowest limit, cal- 

culated using the Optimal Interval Method on the 
spin-independent WIMP-nucleon elastic-scattering cross- 
section for WIMP masses larger than 10 GeV c~^ (see 
Fig. 1381) . These limits begin to constrain some of the 
cosmologically-interesting supersymmetric theories un- 
der the MSSM framework |6l| (see Fig. Of par- 
ticular importance is the exclusion of most of the funnel- 



shaped region in Fig. which corresponds to parame- 
ter space allowed by certain MSSM models where grand 
unification constraints are relaxed |6^ . Other super- 
symmetric frameworks such as mSUGRA '65', '665 ^-^^d 
split-SUSY ,67„ j6Sj remain beyond the sensitivity of this 
WIMP search. These data confirm that events detected 
by CDMS at SUF ^ and those detected by EDEL- 
WEISS 0\ were not a WIMP signal. In addition, these 
data show that the claimed DAMA signal does not 
come from WIMPs of masses greater than 10 GeV c~^ 
distributed in a standard galactic halo and interacting 
through spin-independent interactions. 

In addition to spin-independent (SI) WIMP interac- 
tions, it is also possible that WIMPs may possess spin- 
dependent (SD) interactions with atomic nuclei 0] . Such 
interactions provide additional avenues by which direct 
detection experiments may constrain WIMP models, as 
well as alternative interpretations for reports of positive 
WIMP detection signals jlll . CDMS obtains its SD sen- 
sitivity through the presence of two odd-neutron isotopes 
in its target materials: ^^Si (4.68% of natural Si) and 
^^Ge (7.73% of natural Ge). These nuclides permit us 
to interpret our limits on nuclear-recoil rates in terms of 
limits on WIMP-nucleon spin-dependent cross-sections. 
We proceed under the same assumptions of a standard 
galactic halo, following the method described in 1731 and 
using the form factors given in 't^J for Si and '75] for 
Ge. Figure ITTI shows the resultin g, lim its, calculated us- 
ing the Optimal Interval Method [621 . These data yield 
the lowest WIMP-neutron limit over a substantial range 
of WIMP masses (76,, JT} and are comparable to those 
quoted in [t^, but are ~20% stronger at the minimum 
due to our full knowledge of the detector efficiencies. We 
may also express our limits in terms of allowed regions in 
the Op — an plane of WIMP-nucleon SD couplings, tak- 
ing proper account of the finite-momentum effects de- 
scribed in [ll, [73. The allowed region for WIMPs of 
mass 50 GeV c~^ is shown in Fig. 02 Further details of 
this analysis may be found in [Tq . 



VII. SUMMARY 

From October 11, 2003 until January 11, 2004, the 
CDMS collaboration searched for WIMP dark matter at 
the Soudan Underground Laboratory using four Ge and 
two Si ZIP detectors. The expected background consists 
of 0.7 ± 0.3 betas from contamination on the detector 
surfaces. Out of nearly one million events, only one sur- 
vives all of our analysis cuts. Under the assumptions of a 
standard galactic halo, these data set the world's lowest 
limits on the WIMP-nucleon cross-section in the case of 
spin-independent interactions or spin-dependent interac- 
tions with neutrons. 
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FIG. 35: Events from the ^^^Ci calibrations with recoils between 5-100 keV. Solid lines indicate the 2a nuclear-recoil bands. 
Dashed lines indicate the recoil (vertical) and ionization (curved) analysis thresholds. All the events shown satisfy all of the 
WIMP-search cuts (e.g. single-scatter and fiducial volume) , excluding the ionization-yield and phonon-timing cuts. Gray points 
fail the timing cuts and black points pass the timing cuts. Compare with Fig. 1361 and Fig. 1371 Note that the timing cut rejects 
~20% of the events in the nuclear-recoil band while rejecting >75% of the events in the electron-recoil band. 
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FIG. 36: Unvetoed single-scatters in the fiducial volume for each of the six ZIP detectors of Tower 1. These data correspond 
to the current analysis of 52 livedays of the Run 118 WIMP search at Soudan. Solid lines indicate the 2a nuclear-recoil bands. 
Dashed lines indicate the recoil (vertical) and ionization (curved) analysis thresholds. Gray points fail the timing cuts and 
black points pass the timing cuts. A single event in Z5 with a 64 keV recoil satisfies all of the criteria for a WIMP candidate 
but is consistent with our expected misidentification of electron recoils near the surface. Compare with Fig. Detector Z6 
was not included in this analysis due to a known ^''C contamination. 
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FIG. 37: Unvetoed single-scatters in the fiducial volume for each of the six ZIP detectors of Tower I. These data correspond 
to the initial analysis of 52 livedays of the Run 118 WIMP search at Soudan for which we inadvertently analyzed some of the 
ionization pulses using the time-based fit. Solid lines indicate the 2(t nuclear-recoil-acceptance bands. Dashed lines indicate the 
recoil (vertical) and ionization (curved) analysis thresholds. Gray points fail the timing cuts and black points pass the timing 
cuts. This algorithm led to an overly severe fiducial volume cut which rejected the event in Z5 that passed all of the other cuts. 
Compare with Fig. 1351 Detector Z6 was not included in this analysis due to a known contamination. 
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FIG. 38: Experimental limits in the WIMP-nucleon cross- 
section versus WIMP mass parameter space for spin- 
independent interactions. The cross-section is normalized to 
a single nucleon. The region above the solid (dashed) black 
curves are excluded at 90% confidence by the current (initial) 
analysis of the CDMS-II Soudan WIMP search (Ge detectors). 
The upper thin black curve is the CDMS-II Soudan 90% C.L. 
exclusion limit for the Si detector. The solid gray curve is 
the EDELWEISS i69l| exclusion limit. The dashed gray curve 
is the ZEPLIN I [tI exclusion limit. The 3cr DAMA Fp 
detection region is shown in light gray. Plot courtesy of |7^ . 
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FIG. 39: Comparison of the current (solid curve) and initial 
(dashed curve) CDMS-II exclusion limits with predicted re- 
gions in the WIMP-nucleon cross-section versus WIMP mass 
parameter space for spin-independent interactions. The cross- 
section is normalized to a single nucleon. The light gray re- 
gion corresponds to parameter space that is consistent with 
the standard constraints of the MSSM framework and mea- 
surements of the relic abundance of dark matter ;63| . The two 
dark gray regions correspond to additional allowed parameter 
space under the relaxation of some of the grand unification 
constraints 64]. Plot courtesy of [73| . 



came to fruition. 
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FIG. 40: Comparison of the current (solid curve) and initial FIG. 42: Limits on WIMP-nucleon couplings {up and a„) with 
(dashed curve) CDMS-II exclusion limits with predicted re- a WIMP mass of 50 GeV c . The black curve is the Ge limit 
gions in the WIMP-nucleon cross-section versus WIMP mass ^nd the gray curve is the Si limit, 
parameter space for spin-independent interactions. The cross- 
section is normalized to a single nucleon. The two light 
gray regions correspond to parameter space allowed in the 
mSUGRA framework (65ll66| . The dark gray region and black 
x's correspond to allowed regions under the assumptions of 
split-SUSY miH. Plot courtesy of [72. 
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FIG. 41: CDMS-II (Tower 1) spin-dependent limits. Black 
lines are Ge limits, gray lines are Si limits. Dotted lines are 
proton limits and solid lines are neutron limits. 
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